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ABSTRACT 


Due  the  remarkable  growth  in  the  size  and  complexity  of  airlift  operations,  there  is  an 
increased  need  for  planning  tools  to  assist  decision  makers  with  issues  ranging  from  selecting 
the  number  and  types  of  aircraft  for  an  airlift  fleet  to  making  informed  decisions  with  respect 
to  investing  or  divesting  in  overseas  air  bases.  In  Fiscal  Year  (FY)  94  research  was  initiated 
in  the  Operations  Research  Department  of  the  Naval  Postgraduate  School  in  response  to  a 
request  from  the  United  States  Air  Force  Studies  and  Analyses  Agency  and  resulted  in  the 
development  of  a  high  fidelity  strategic  airlift  optimization  model  called  Throughput  II.  The 
model  is  formulated  as  a  multi-period,  multi-commodity  linear  programming  model  for 
determining  the  maximum  on-time  throughput  of  cargo  and  passengers  that  can  be 
transported  with  a  given  fleet  or  given  network,  subject  to  appropriate  physical  and  policy 
constraints.  Troop  and  equipment  movement  requirements  are  specified  by  the  Time  Phase 
Force  Deployment  Data  (TPFDD).  An  optimization  model  that  utilizes  the  full  level  of  detail 
available  in  a  TPFDD  would  be  of  intractable  size.  Moreover,  it  is  not  necessary  to  build  a 
model  with  such  a  fine  level  of  detail  in  order  to  obtain  the  important  insights  required  to 
assist  decision  makers.  Therefore  Throughput  II  replaces  the  potentially  large  set  of  airfields 
with  a  smaller  set  of  centroids  and  schedules  aircraft  through  these  aggregated  airfields. 
Currently  route  selection  is  performed  manually,  by  an  expert,  who  incorporates  a  variety  of 
factors  based  on  his/her  experience.  In  this  thesis  we  develop  techniques  for  selecting  a  set 
of  candidate  routes  for  any  deployment  scenario  without  requiring  historical  data  or  extensive 
interaction  with  an  expert.  An  analyst  should  be  concerned  about  two  potentially  detrimental 
effects  of  these  preprocessing  procedures.  First,  infeasibility  may  be  introduced  by 
aggregation  and  second,  Throughput  II  may  provide  suboptimal  solutions  since  we  consider 
a  limited  number  of  routes.  To  address  these  issues,  a  postprocessing  step  can  be  used  to 
screen  for  constraint  violations  and  to  perform  sensitivity  analysis  with  respect  to  alternative 
routing  options. 
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EXECUTIVE  SUMMARY 


The  National  Defense  Authorization  for  Fiscal  Year  (FY)  1991  required  the 
Department  of  Defense  to  conduct  a  comprehensive  study  to  determine  the  mobility 
requirements  for  the  United  States  Armed  Forces,  The  Mobility  Requirements  Study 
(MRS)  examined  all  aspects  of  the  mobility  question,  from  domestic  transportation  to 
intertheater  to  intratheater  requirements.  The  MRS  analyzes  sea,  air  and  amphibious  lift, 
surface  transportation  and  prepositioning  requirements  to  provide  Congress  with  an 
integrated  plan  for  procuring  the  necessary  lift  assets  and  infrastructure  for  power 
projection  in  the  21st  century.  In  January  1991,  the  Force  Structure,  Resource  and 
Assessment  Directorate  (J8)  of  the  Joint  Staff  assumed  responsibility  for  the  MRS.  The 
Mobility  Optimization  Model  (MOM)  [Wing  et  al,  1991]  was  developed  by  J8  with 
assistance  from  the  Naval  Postgraduate  School  (NPS)  within  a  short  period  of  six  weeks. 
MOM  was  developed  to  determine  airlift,  sealift  and  prepositioning  requirements.  It 
modeled  single  onload  and  offload  location  with  a  time  window  for  delivery.  Even  though 
MOM  produces  reasonable  high  level  results,  it  is  not  suitable  for  answering  the  mobility 
questions  currently  sought  by  USAF  Studies  and  Analyses  Agency  (USAF/SAA)  because 
it  does  not  model  airlift  assets,  infrastructure  and  the  associated  constraints  in  sufficient 
detail. 

A  second  linear  programming  (LP)  model  for  strategic  airlift  was  USAF/SAA’s 
Throughput  Model  (Yost,  1 994).  This  model  was  developed  to  determine  an  optimal 


airlift  fleet  mix  based  on  given  movement  requirements  and  to  compare  different  airlift 
fleets  in  terms  of  cargo  throughput  capabilities.  Throughput  proved  to  be  a  useful  model 
but  had  some  inherent  limitations  due  to  its  static  nature  with  respect  to  time  Important 
system  details  such  as  delivery  time  windows  for  specific  units  and  unit  closure  dates  are 
not  incorporated  In  FY  94  research  was  performed  at  NFS  in  response  to  a  request  from 
the  USAF/SAA,  and  resulted  in  the  development  of  a  higher  fidelity  strategic  airlift  model 
called  Throughput  II  (Lim,  1994)  which  has  since  been  enhanced  (Morton,  et  a/.,  1995).  In 
this  research,  the  strategic  airlift  assets  optimization  problem  is  formulated  as  a  multi¬ 
period,  multi-commodity  network-based  linear  programming  model,  with  a  large  number 
of  side  constraints.  The  objective  of  the  model  is  to  minimize  late  deliveries  subject  to 
physical  and  policy  constraints,  such  as  aircraft  utilization  limits  and  airfield  handling 
capacities. 

A  typical  Time  Phased  Force  Deployment  Data  (TPFDD)  document  specifies  over 
14000  specific  movement  requirements  over  180  days  from  120  Aerial  Ports  of 
Embarkation  (APOE)  to  90  Aerial  Ports  of  Debarkation  (APOD)  through  a  network 
infrastructure  with  a  large  number  of  potential  enroute  stops.  To  attempt  to  build  an 
optimization  model  that  incorporates  the  full  level  of  detail  available  in  a  TPFDD  would 
result  in  a  mathematical  model  of  intractable  size.  Moreover,  it  is  not  necessary  to  build  a 
model  with  such  a  fine  level  of  detail  in  order  to  obtain  the  important  insights  required  to 
assist  decision  makers.  Therefore,  Throughput  and  Throughput  II  both  use  data 
aggregated  from  the  TPFDD  for  input.  Both  models,  for  example,  replace  the  potentially 


large  set  of  airfields  with  a  smaller  set  of  centroid  airfields  and  schedule  aircraft  through 
these  aggregate  airfields. 

We  develop  a  mixed  integer  location-allocation  centroid  model  that  may  be  used  as 
a  preprocessing  step  Aggregation  of  airfields  may  lead  to  infeasibilities  in  the  resulting 
model  with  respect  to  ability  to  fly  certain  flight  legs  due  to  altered  inter-airfield  distances. 
This  is  the  principle  reason  that  we  propose  a  centroid  aggregation  model  that  has  the 
primary  objective  of  minimizing  the  maximum  centroid-satellite  distance.  A  pure  minimax 
model,  will  typically  have  many  multiple  optimal  solutions.  Some  optima,  however,  are 
better  than  others  with  respect  to  criteria  ignored  by  the  minimax  objective  (e.g.,  the 
satellite-centroid  for  second  furthest  airfield).  Because  of  this  drawback,  a  secondary 
minisum  objective  term  is  included  in  the  centroid  selection  model. 

Currently,  route  selection  is  performed  manually,  by  an  expert,  who  incorporates  a 
variety  of  factors  including  critical  legs  for  specific  aircraft  types,  aircrew  flight  time 
restrictions  and  commercial  track  avoidance  of  certain  countries.  Such  a  procedure  is  time 
consuming  and  may  be  adequate  for  flying  to  destinations  with  which  USAF  has  extensive 
experience.  However,  more  difficulties  may  arise  when  flying  to  an  unfamiliar  destination. 
We  develop  heuristic  techniques  for  selecting  the  optimal  routes  for  any  deployment 
scenario  without  requiring  historical  data  and  extensive  interaction  with  an  expert.  The 
route  selection  module  is  implemented  in  the  Pascal  programming  language.  Feasible 
routes  are  sorted  with  respect  to  their  lengths  for  each  origin-destination  pair.  Inefficient 
and  infeasible  routes,  due  to  constraints  such  as  crew  flying  hours  or  crew  stages,  are  then 
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eliminated.  This  automated  process  has  the  advantage  that  routes  may  be  generated 
without  requiring  the  extensive  experience  and  knowledge  of  an  expert.  Nevertheless, 
there  are  certain  system  requirements  (e  g.,  country  permissions)  that  are  difficult  to 
embed  in  such  a  program,  and  expert  knowledge  may  still  be  required  to  rule  out  the 
feasibility  of  certain  routes. 

The  elimination  and  preprocessing  step  may  still  generate  a  large  number  of  routes, 
and  it  may  not  be  computationally  practical  to  include  all  the  candidate  routes  in  the 
Throughput  II  model  For  this  reason,  a  method  is  applied  to  convert  the  large  problem 
into  a  smaller  problem  of  manageable  size.  This  reduction  is  achieved  by  solving  the  linear 
program  over  a  subset  of  its  columns;  the  result  of  which  may  be  suboptimal.  Through 
standard  LP  sensitivity  analysis  we  can  test  whether  a  certain  route,  not  in  the  original 
subset,  is  attractive.  If  it  is  attractive,  it  can  be  put  in  the  basis  as  a  new  column  and  is  used 
to  improve  the  route  generation  model.  This  procedure  can  be  automated  and  is  known  as 
column  generation;  in  our  case,  the  route  selection  model  will  be  the  column  generator. 
Including  new  attractive  routes  and  iterating  the  Throughput  II  model  can  be  continued 
until  a  satisfactory  solution  is  achieved. 
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INTRODUCTION 


In  spite  of  the  sophisticated  inventories  and  impressive  capabilities  of  modern  air 
forces,  military  aviation  is  still  the  youngest  of  all  major  military  operations.  Except  for  a 
few  early  attempts,  such  as  using  balloons  and  airships  as  instruments  of  war,  air  power 
has  been  an  extremely  effective  means  of  achieving  military  objectives,  mainly  because  of 
its  great  versatility,  for  the  last  80  years.  Among  their  many  missions,  airforces  provide 
two  key  tasks:  delivering  firepower  directly  by  air-to-surface  and  air-to-air  weapons,  and 
rapidly  deploying  combat  units  over  substantial  distances.  While  the  use  of  military  air 
power  is  relatively  new,  air  mobility  operations  are  of  even  more  recent  vintage.  The 
substantial  exploitation  of  military  airlift  dates  only  from  1938.  Prior  to  that  time,  except 
for  some  primitive  transportation  efforts  on  a  limited  scale  via  air  transport  systems,  the 
deployment  and  logistic  support  of  military  forces  was  undertaken  mostly  via  surface 
transport  systems. 

There  are  three  fiindemental  disadvantages  of  surface  deployment,  which  cannot  be 
eliminated  by  speed  and  efficiency  provided  by  the  newly  developing  ideas  and  technology 
of  this  century.  Firstly,  surface  deployment  over  substantial  distances  is  slower  than  air 
transport  systems.  Secondly,  surface  deployment  is  heavily  restricted  by  geographical 
constraints.  The  third,  and  most  important  reason,  concerns  its  susceptibility  to  enemy 
attacks.  On  the  other  hand,  the  airlift  has  less  comparative  capacity  than  the  sealift  does. 
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but  because  of  its  speed  and  range,  all  the  modern  armed  forces  have  incorporated  airlift 
capability  into  their  transportation  systems. 

During  World  War  II,  in  addition  to  the  effectiveness  of  fighting,  attacking  and 
bombing  units  of  air  forces,  air  transport  operations  proved  to  be  an  efficient  means  of 
projecting  forces  to  strategic  locations.  This  added  another  aspect  to  the  obvious  potential 
of  air  forces  as  a  direct  force,  namely  the  use  of  air  forces  as  an  indirect  force.  The 
importance  of  this  latter  feature  of  air  forces  is  easily  recognized  by  the  success  of  airlift 
operations  which  were  totally  impossible  by  other  means  in  the  same  range  and  time 
frames. 

Despite  prepositioned  fuel,  ammunition,  and  equipment,  the  magnitude  of  the 
United  Nations/  United  States  airlift  effort  during  Desert  Shield  and  Desert  Storm  was 
unprecedented  By  10  March  1991,  strategic  airlift  had  moved  over  500,000  people  and 
540,000  tons  of  cargo  to  the  theater.  At  the  height  of  the  Desert  Shield  airlift.  Military 
Airlift  Command’s  (MAC’s)  cargo  movement  averaged  17  million  ton-miles  per  day’.  By 
comparison,  during  the  1973  Arab-Israeli  War,  US  airlift  moved  4.4  million  ton-miles  per 
day  Other  historical  comparisons  include  the  World  War  II  “Hump”  at  0.9  million  ton- 
miles  per  day,  the  Berlin  Airlift  at  1 .7  million  ton-miles  per  day,  and  Operation  Just  Cause 
at  2.0  million  ton-miles  per  day  (Gulf  War  Air  Power  Survey,  1993). 


'  US  strategic  airlift  requirements  and  capabilities  may  be  measured  in  terms  of  million-ton-miles  per 
day.  i  c.,  in  multiples  of  the  capacity  to  move  one  ton  of  cargo  by  air  a  distance  of  one  mile  in  one  day. 
This  standard  of  measurement,  though  a  useful  planning  tool  and  gauge  of  present  and  future  strategic 
airlift  capabilities,  does  not,  of  course,  take  into  account  such  real  world  constraints  as  aircraft  and  crew 
exhaustion,  weather,  availability  of  airfields,  overflight  rights,  and  possible  enemy  action. 
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Due  to  the  remarkable  growth  in  the  size  and  complexity  of  airlift  operations,  there 
is  increased  need  for  planning  tools  to  assist  decision  makers  with  issues  ranging  from 
selecting  the  number  and  types  of  aircraft  for  an  airlift  fleet  to  making  informed  decisions 
with  respect  to  investing  or  divesting  in  overseas  air  bases.  Although  it  must  be  interpreted 
carefully,  this  “search  for  optimality”  is  a  very  important  theme  in  Operations  Research. 
The  field  of  Operations  Research  provides  tools  for  attempting  to  find  the  best,  or  optimal, 
solution  to  mathematical  models  of  the  problem  under  consideration. 

A.  BACKGROUND 

The  National  Defense  Authorization  for  Fiscal  Year  (FY)  1991  required  the 
Department  of  Defense  to  conduct  a  comprehensive  study  to  determine  the  mobility 
requirements  for  the  United  States  Armed  Forces.  The  Mobility  Requirements  Study 
(MRS)  examined  all  aspects  of  the  mobility  question,  from  domestic  transportation  to 
intertheater  to  intratheater  requirements.  The  MRS  analyzes  sea,  air  and  amphibious  lift, 
surface  transportation  and  prepositioning  requirements  to  provide  Congress  with  an 
integrated  plan  for  procuring  the  necessary  lift  assets  and  infrastructure  for  force 
projection  in  the  21  st  century.  In  January  1 991,  the  Force  Structure,  Resource  and 
Assessment  Directorate  (J8)  of  the  Joint  Staff  assumed  responsibility  for  the  MRS.  The 
Mobility  Optimization  Model  (MOM)  (Wing  et  al.,  1991)  was  developed  by  J8  with 
assistance  from  the  Naval  Postgraduate  School  (NPS)  within  a  short  period  of  six  weeks. 
MOM  was  developed  to  determine  airlift,  sealift  and  prepositioning  requirements.  For 
airlift,  it  modeled  single  onload  and  offload  locations  with  a  time  window  for  delivery. 
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Even  though  MOM  produces  reasonable  high  level  results  it  is  not  suitable  for  answering 
the  mobility  questions  currently  sought  by  USAF  Studies  and  Analyses  Agency 
(USAF/SAA)  because  it  does  not  model  airlift  assets,  infrastructure  and  the  associated 
constraints  in  sufficient  detail. 

A  second  linear  programming  (LP)  model  for  strategic  airlift  was  USAF/SAS’s 
Throughput  Model  (Yost,  1994).  This  model  was  developed  to  determine  an  optimal 
airlift  fleet  mix  based  on  given  movement  requirements  and  to  compare  different  airlift 
fleets  in  terms  of  cargo  throughput  capabilities.  Throughput  proved  to  be  a  useful  model 
but  had  some  inherent  limitations  due  to  its  static  nature  with  respect  to  time.  Important 
system  details  such  as  delivery  time  windows  for  specific  units  and  unit  closure  dates  are 
not  incorporated.  In  FY  94,  research  was  performed  at  NPS  in  response  to  a  request  from 
the  USAF/SAA,  and  culminated  with  development  of  a  higher  fidelity  strategic  airlift 
model  called  Throughput  II  (Lim,  1994).  In  this  research,  the  strategic  airlift  assets 
optimization  problem  is  formulated  as  a  multi-period,  multi-commodity  network -based 
linear  programming  model,  with  a  large  number  of  side  constraints.  The  objective  of  the 
model  is  to  minimize  late  deliveries  subject  to  physical  and  policy  constraints,  such  as 
aircraft  utilization  limits  and  airfield  handling  capacities. 

B.  PROBLEM  STATEMENT 

Executing  a  deployment  -  delivering  units  and  equipment  to  theater  efficiently  and 
effectively  -  with  limited  aircraft  and  infrastructure  assets  is  a  formidable  task.  The 
National  Command  Authority  (NCA),  the  Joint  Chief  of  Staff,  and  major  field 
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commanders  have  access  to  the  Worldwide  Military  Command  and  Control  System 
(WWMCCS),  which  includes  a  command  and  control  system  called  Joint  Operational 
Planning  and  Execution  System  (JOPES).  The  Operational  Plan  (OPLAN)  is  entered  as  a 
requirement  and  JOPES  generates  a  Time  Phased  Force  Deployment  List  (TPFDL).  This 
list  called  the  Time  Phased  Force  Deployment  Data  (TPFDD)  document,  and  contains 
information  pertaining  to  the  force  and  movement  requirements  of  the  OPLAN.  With  this 
information,  units  able  to  supply  force  requirement  can  be  identified  and  transportation 
arrangements  can  be  made.  A  typical  TPFDD  specifies  over  14000  specific  movement 
requirements  over  180  days  from  120  Aerial  Ports  of  Embarkation  (APOE)  to  90  Aerial 
Ports  of  Debarkation  (APOD)  through  a  network  infrastructure  with  a  large  number  of 
potential  enroute  stops  To  attempt  to  build  an  optimization  model  that  utilizes  the  full 
level  of  detail  available  in  a  TPFDD  would  result  in  a  mathematical  model  of  intractable 
size  Moreover,  it  is  not  necessary  to  build  a  model  with  such  a  fine  level  of  detail  in  order 
to  obtain  the  important  insights  required  to  assist  decision  makers.  Therefore,  Throughput 
and  Throughput  II  both  use  data  aggregated  from  the  TPFDD  for  input.  Both  models,  for 
example,  replace  the  potentially  large  set  of  airfields  with  a  smaller  set  of  centroid  airfields 
and  schedule  aircraft  through  these  aggregate  airfields. 

A  mixed  integer  programming  model  for  centroid  selection  has  been  developed  by 
AFSAA  that  may  be  used  as  a  preprocessing  step  for  either  model  (Yost,  1994).  The  basic 
idea  in  this  model  is  to  first  define  candidate  centroid  airfields  based  on  their  passenger  and 
cargo  flows  and  then  to  solve  the  model  that  minimizes  a  weighted  sum  of  distances 
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between  the  centroids  and  their  satellite  airfields.  A  primary  shortcoming  of  such  an 
aggregation  procedure  is  that  it  may  lead  to  infeasible  solutions  in  the  subsequent  airlift 
optimization  model.  Airfields  aggregated  into  centroids  may  inappropriately  “borrow” 
capacity  from  one  another.  As  a  result,  it  may  be  discovered  on  disaggregation  that  the 
solution  to  Throughput  II  has  assigned  too  many  sorties  to  an  airfield.  Such  an  outcome 
may  occur,  for  example,  when  a  relatively  small  capacity  airfield  with  a  large  flow 
requirement  is  aggregated  with  another  airfield  that  has  a  large  capacity  value.  A  second 
potential  drawback  in  the  current  method  of  aggregating  airfields  is  that  the  centroids  are 
selected  to  minimize  the  sum  of  weighted  distances  to  all  corresponding  satellite  airfields, 
and  the  distances  between  airfields  in  Throughput  II  are  determined  from  the  centroid  not 
the  satellite.  As  a  result,  there  are  cases  in  which  a  particular  aircraft  type  may  be  capable 
of  flying  a  route  which  originates  from  the  centroid  but  not  from  the  satellite  airfield. 
Conversely,  certain  routes  may  be  eliminated  unnecessarily  due  to  airfield  aggregation.  As 
a  result  of  these  aggregation  effects,  the  fidelity  of  the  associated  Throughput  II  model 
may  be  reduced. 

Currently,  route  selection  is  performed  manually,  by  an  expert,  who  incorporates  a 
variety  of  factors  including  minimum  payload  ferry  ranges^  for  specific  aircraft  types, 
aircrew  flight  time  restrictions  and  commercial  track  avoidance  of  certain  countries.  Such 
a  procedure  is  time  consuming  but  may  be  adequate  for  flying  to  destinations  with  which 


‘  This  is  the  longest  distance  which  an  aircraft  can  fly  with  minimum  reasonable  payload  (20  %  of  its 
maximum  capacity)  without  refueling. 


6 


USAF  has  extensive  experience;  however,  difficulties  may  arise  when  flying  to  an 
unfamiliar  destination. 

In  summary,  the  size  of  the  infrastructure  network  is  currently  reduced  by  a 
preprocessing  step  that  selects  airfields  via  the  centroid  model  and  eliminates,  by  hand,  a 
large  number  of  infeasible  or  inefficient  routes  The  results  of  Throughput  II  have 
indicated  that  for  some  important  deployment  scenarios  the  airlift  system  may  be 
infrastructure  constrained,  in  part,  due  to  the  aircraft  handling  capability  of  certain  enroute 
airfields.  This  makes  the  process  of  preselecting  centroid  airfields  as  well  as  delivery  and 
recovery  routes  very  important.  It  is  desirable  to  minimize  the  effect  of  capacity 
borrowing  Moreover,  it  may  be  possible  to  improve  the  system’s  throughput  by  proper 
aircraft-route  scheduling.  An  automated  system  that  can  generate  routes  for  given  origin- 
destination  pairs  could  dramatically  aid  the  analyst  in  the  selection  of  deployment  and 
recovery  routes. 

C.  METHODOLOGY 

This  thesis  will  consist  of  two  main  parts  as  illustrated  in  Figure  1.1.  The  first  part 
is  a  preprocessing  procedure  that  selects  a  set  of  centroids  and  candidate  routes  that  will 
be  used  as  input  to  Throughput  II  or  any  other  appropriate  optimization  model.  The 
second  part  is  a  postprocessing  procedure  that  seeks  to  determine  whether  the 
Throughput  II  solution  may  be  disaggregated  into  a  valid  solution  for  the  original  model 
and/or  improved  by  providing  Throughput  II  with  additional  routing  options. 
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Figure  1.1.  Route  Selection  Model 


1.  Preprocessing  Modules 

The  preprocessing  step  consists  of  two  distinct  modules;  first  a  model  is  solved  to 
determine  the  centroids  and  then  a  heuristic^  generates  routes  from  the  origin  centroid 
airfields  through  enroute  airfields,  to  the  destination  centroid  airfields  for  all  relevant 
aircraft  types. 


CL  Centroid  Model 

We  develop  a  location-allocation  model  with  a  primary  objective  to 
minimize  the  maximum  centroid-satellite  distance  and  a  secondary  objective  to  minimize 
the  sum  of  the  corresponding  distances.  Aggregation  of  airfields  may  lead  to  infeasibilities 


^  “A  heuristic  is  a  rule  of  thumb,  strategy,  or  trick  used  to  improve  the  efficiency  of  a  system  which  tries  to 
discover  the  solutions  of  complex  problems.”(Slagle,  1971). 
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in  the  resulting  model,  with  respect  to  an  aircraft’s  ability  to  fly  certain  flight  legs,  due  to 
altered  inter-airfield  distances.  This  is  the  primary  reason  that  we  propose  a  centroid 
aggregation  model  that  has  a  primary  objective  of  minimizing  the  maximum  centroid- 
satellite  distance  A  pure  minimax  model,  typically  has  many  multiple  optimal  solutions. 
Some  optima,  however,  are  better  than  others  with  respect  to  criteria  ignored  by  the 
minimax  objective  (e  g  ,  the  satellite-centroid  for  second  furthest  airfield).  Because  of  this 
drawback,  the  secondary  minisum  objective  term  is  included  in  the  centroid  selection 
model 

h.  Route  Selection  Model 

The  route  selection  module  is  implemented  in  the  Pascal  programming 
language  Feasible  routes  are  sorted  with  respect  to  their  lengths  for  each  origin- 
destination  pair.  Inefficient  and  infeasible  routes,  due  to  constraints  such  as  crew  flying 
hours  or  crew  stages,  are  then  eliminated.  This  automated  process  has  the  advantage  that 
routes  may  be  generated  without  requiring  the  extensive  experience  and  knowledge  of  an 
expert.  Nevertheless,  there  are  certain  system  requirements  (e.g.,  country  permissions) 
that  are  difficult  to  embed  in  such  a  program,  and  expert  knowledge  may  still  be  required 
to  rule  out  the  feasibility  of  certain  routes. 

2.  Postprocessing  Module 

The  final  output  of  the  preprocessing  step  is  a  set  of  candidate  routes  that  can  be 
used  as  input  data  for  Throughput  II  or  any  other  appropriate  model  such  as  the 
Stochastic  Airlift  Optimization  Model  (Goggins,  1995).  However  the  preprocessing  step 
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may  still  generate  a  large  number  of  routes,  and  it  may  not  be  computationally  practical  to 
include  all  the  candidate  routes  in  the  Throughput  II  model.  For  this  reason,  a  method  is 
applied  to  convert  this  large  problem  into  a  smaller  problem  of  manageable  size.  This 
reduction  is  achieved  by  solving  the  linear  program  over  a  subset  of  its  columns,  the  result 
may  be  suboptimal.  Through  standard  LP  sensitivity  analysis  we  can  test  whether  a  certain 
route,  not  in  the  original  subset,  is  attractive.  If  it  is  attractive,  it  can  be  put  in  the  basis  as 
a  new  column  and  is  used  to  improve  the  route  generation  model.  This  procedure  can  be 
automated  and  is  known  as  column  generation,  in  our  case,  the  route  selection  model  will 
be  the  column  generator.  This  process  of  including  new  attractive  routes  in  the 
Throughput  II  model,  solving  the  corresponding  LP,  and  testing  for  optimality  can  be 
continued  until  a  solution  to  the  original  problem  is  achieved. 

D.  SUMMARY 

In  this  thesis  we  develop  techniques  for  selecting  the  optimal  routes  for  any 
deployment  scenario  without  requiring  historical  data  and  extensive  interaction  with  an 
expert.  The  results  of  Throughput  II  have  indicated  that  for  some  important  deployment 
scenarios  the  airlift  system  may  be  infrastructure  constrained.  This  makes  the  process  of 
preselecting  centroid  airfields  very  important.  In  the  first  part  of  this  thesis  we  develop  an 
efficient  procedure  for  solving  the  combination  minimax-minisum  facility  location  model 
to  determine  centroid  airfields  that  will  help  curb  detrimental  aggregation  effects. 
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An  analyst  should  be  concerned  about  two  potentially  detrimental  effects  of  the 
preprocessing  procedure  Either  infeasibility  may  be  introduced  by  aggregation  or 
Throughput  II  may  yield  a  suboptimal  solution  since  only  a  limited  number  of  routes  are 
considered  To  address  these  issues,  a  postprocessing  step  includes 
i.  disaggregation  and  feasibility  tests, 

ii  sensitivity  analysis  with  respect  to  alternate  routing  options. 

Together,  these  measures  assist  Throughput  II  in  providing  tractible,  as  well  as 
valid  results. 
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11. 


REVIEW  OF  THROUGHPUT  11  MODEL 


The  Throughput  II  strategic  airlift  optimization  model  was  developed  by  Captain 
Lim  Teo  Weng  (Lim,  1994)  and  has  since  been  enhanced  in  ongoing  research  at  the  Naval 
Postgraduate  School  (Morton  et  al,  1995).  The  mathematical  formulation  of  Throughput 
II  is  contained  in  Appendix  A.  The  Throughput  II  model  is  a  multi-period  model  for 
determining  the  maximum  on-time  throughput  of  cargo  and  passengers  that  can  be 
transported  with  a  given  fleet  over  a  given  network,  subject  to  appropriate  physical  and 
policy  constraints  The  objective  function  minimizes  the  total  weighted  penalties  for  late 
deliveries  and  nondeliveries. 

The  preprocessing  and  postprocessing  techniques  developed  in  this  thesis  are 
utilized  by  the  Throughput  II  optimization  model  as  indicated  in  Chapter  I.  For  this 
reason,  we  provide  a  review  of  Throughput  II  in  this  chapter. 

A.  MODEL  FEATURES 

The  Throughput  II  model  has  been  designed  to  handle  many  of  the  airlift  system’s 
particular  features  and  modes  of  operation.  The  model  is  a  strategic  airlift  model,  meaning 
that  it  considers  inter-theater  but  not  intra-theater  deliveries.  The  major  features  of  the 
airlift  system  currently  captured  by  the  model  include  (Morton  etal,  1995): 

•  Multiple  origins  and  destinations:  The  model  routes  aircraft  through  multiple  origin, 
enroute  and  destination  airfields. 
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Flexible  routing  structure:  The  air  route  structure  supported  by  the  model  includes 
delivery  and  recovery  routes  with  a  variable  number  of  enroute  stops  (usually  between 
zero  and  three)  This  gives  the  model  the  option  of  short-range  flights  with  heavier 
loads  or  long-range  flights  with  lighter  loads.  For  further  routing  flexibility,  the  model 
also  allows  the  same  aircraft  to  fly  different  delivery  and  recovery  routes  on  the  same 
mission 

Aircraft-to-route  restrictions:  The  user  may  impose  aircraft-to-route  restrictions,  e  g., 
military  aircraft  may  only  use  military  airfields  for  enroute  stops.  This  particular 
provision  arises  because  the  USAF  Air  Mobility  Command  (AMC)  may  call  upon 
civilian  commercial  airliners  to  augment  USAF  aircraft  in  a  deployment,  under  the 
Civil  Reserve  Airlift  Fleet  (CRAF)  program.  The  model  distinguishes  between  USAF 
and  CRAF  aircraft. 

Aircraft  assets  can  be  added  over  time.  This  adds  realism  to  the  model  because  CRAF 
and  other  aircraft  may  take  time  to  mobilize  and  are  typically  unavailable  at  the  start  of 
a  deployment. 

Delivery  time  windows:  In  a  deployment,  a  unit  is  ready  to  move  on  its  available-to- 
load  date  (ALD)  and  has  to  arrive  at  the  theater  by  its  required-delivery-date  (RDD). 
This  aspect  of  the  problem  has  been  incorporated  in  the  model  through  user-specified 
time  windows  for  each  unit.  The  model  treats  the  time  window  as  “elastic”  in  that 
cargo  may  be  delivered  late,  subject  to  a  penalty. 


ASSUMPTIONS 


The  assumptions  used  in  the  model  are  as  follows: 

Airfield  capacity  is  represented  by  a  single  aggregate  figure,  called  Maximum-on  - 
Ground  (MOG).  The  literal  translation  of  MOG  as  the  maximum  number  of  planes 
that  can  be  simultaneously  on  the  ground  at  an  airfield  is  misleading,  because  MOG  is 
used  to  convey  more  than  just  the  number  of  parking  spaces.  In  actuality,  airfield 
capacity  depends  on  many  additional  factors  such  as  availability  of  material  handling 
equipment  and  various  ground  services  capacities.  Unfortunately,  data  are  not 
currently  available  to  support  a  multidimensional  MOG  modeling  enhancement. 

Inventoried  aircraft  at  origin  and  destination  airfields  are  considered  not  to  affect  the 
aircraft  handling  capacity  of  the  airfield.  This  assumption  is  not  strictly  valid  since  an 
inventoried  aircraft  takes  up  parking  space  even  if  it  is  not  consuming  services. 


•  Deterministic  ground  time;  Aircraft  turnaround  times  for  onloading  and  offloading 
cargo  and  enroute  refueling  are  assumed  to  be  known  constants,  although  they  are 
naturally  stochastic  This  ignores  the  fact  that  deviations  from  the  given  service  time 
can  cause  congestion  on  the  ground  To  offset  the  optimism  of  this  assumption,  an 
efficiency  factor  is  used  in  the  formulation  of  airfield  capacity  constraints. 

C.  CONCEPTUAL  MODEL  FORMULATION 

The  primary  decision  variables  are  the  number  of  sorties  initiated,  and  the  amount 
of  cargo  and  passengers  carried,  for  each  unit,  by  each  aircraft  type,  via  each  available 
route,  in  each  time  period.  Additional  variables  are  defined  for  the  recovery  flights,  for 
aircraft  inventoried  at  airfields,  and  for  the  possibility  (at  high  penalty  cost)  of  not 
delivering  required  cargo  or  passengers. 

1.  Objective  Function 

The  objective  function  minimizes  the  total  weighted  penalties  for  late  deliveries 
and  nondeliveries  subject  to  appropriate  physical  and  policy  constraints.  The  penalties  are 
weighted  according  to  two  factors:  the  priority  of  the  unit  and  the  degree  of  lateness.  The 
penalty  increases  with  the  amount  of  time  late,  and  non-delivery  has  the  most  austere 
penalty. 

The  anticipated  use  of  the  model  is  for  situations  when  the  given  airlift  resources 
are  insufficient  for  making  all  the  required  deliveries  on  time.  On  the  other  hand,  if  there 
are  enough  resources  for  complete  on-time  delivery,  then  the  model’s  secondary  objective 
function  is  to  choose  a  feasible  solution  that  maximizes  unused  aircraft.  The  motivation  for 
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the  secondary  objective  is  that  if  the  available  aircraft  are  used  as  frugally  as  possible, 
while  still  meeting  the  known  demands  and  observing  the  known  constraints,  then  the 
mobility  system  will  be  as  well  prepared  as  possible  for  unplanned  breakdowns  and 
unforeseen  requirements,  such  as  an  additional  nearly  simultaneous  regional  contingency. 

2.  Constraints 

The  model’s  constraints  can  be  grouped  into  five  categories:  demand  satisfaction, 
aircraft  balance,  aircraft  capacity,  aircraft  utilization,  and  airfield  handling  capacity. 


•  Demand  Satisfaction  Constraints:  The  cargo  demand  constraints  attempt  to  ensure  for 
each  unit  that  the  correct  amount  of  cargo  moves  to  the  required  destination  within  the 
specified  time  window.  The  passenger  demand  constraints  do  the  same  for  each  unit’s 
personnel.  The  demand  constraints  have  elastic  variables  for  late  delivery  and  non¬ 
delivery  The  optimization  will  seek  to  avoid  these  options  if  it  is  possible  with  the 
available  assets,  or  to  minimize  them  if  not. 

•  Aircraft  Balance  Constraints:  These  constraints  keep  physical  count  of  aircraft  by  type 
(e  g.,  C17,  C5,  C141,  etc.)  in  each  time  period.  They  ensure  that  the  aircraft  assets  are 
used  only  when  they  are  available. 

•  Aircraft  Capacity  Constraints:  There  are  three  different  kinds  of  constraints  on  the 
physical  limitations  of  aircraft  —  troop  carriage  capacity,  maximum  payload,  and  cabin 
floor  space  —  which  must  be  observed  at  all  times. 

•  Aircraft  Utilization  Constraints:  These  constraints  ensure  that  the  average  flying  hours 
consumed  per  aircraft  per  day  are  within  AMC’s  established  utilization  rates  for  each 
aircraft  type. 

•  Aircraft  Handling  Capacity  at  Airfields:  These  constraints  ensure  that  the  number  of 
aircraft  routed  through  each  airfield  each  day  is  within  the  airfield’s  handling  capacity. 
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D.  LIMITATIONS 

Two  limitations  of  the  optimization  model  are  the  inability  to  handle  local  traffic 
congestion  and  low  fidelity  with  respect  to  certain  aspects  of  the  airlift  system. 

1.  Time  Resolution 

With  one-day  time  periods,  the  model  can  route  aircraft  in  a  manner  that  causes 
local  congestion.  Another  limitation  is  that  the  model  rounds  the  time  (to  the  nearest  day) 
at  which  an  aircraft  arrives  at  enroute  and  destination  airfields. 

2.  Aircraft  Reliability 

Aircraft  in  need  of  repair  can  have  an  immediate  impact  on  throughput  capability, 
especially  when  airfield  capacity  resources  are  limited. 

3.  Deterministic  Ground  Time 

The  aircraft  ground  times  used  in  the  calculation  of  MOG  consumption  represent 
the  expected  times  for  onload/refuel/offload,  resulting  in  optimistic  throughput  capability. 

4.  Airfield  Aggregation 

Throughput  II  uses  data  aggregated  from  the  TPFDD  for  input.  It  replaces  the 
large  set  of  airfields  with  a  smaller  set  of  centroid  airfields  and  schedule  aircraft  through 
these  aggregate  airfields. 
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5. 


Limited  Route  Structure 


Route  selection  is  performed  manually,  by  an  expert,  who  incorprates  a  variety  of 
factors  including  minimum  payload  ferry  ranges  for  specific  aircraft  types,  aircrew  flight 
time  restrictions  and  commercial  track  avoidance  of  certain  countries.  Such  a  procedure  is 
time  consuming  but  may  be  adequate  for  flying  to  destinations  with  which  USAF  has 
extensive  experience;  however  difficulties  may  arise  when  flying  to  an  unfamiliar 
destination 

Airfield  aggregation  and  limited  route  structure  are  examined  in  this  thesis. 
Limitations  involving  aircraft  reliability  and  stochastic  ground  times  are  examined  by 
Goggins  (1995). 
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III. 


PREPROCESSING  MODULES 


The  preprocessing  step  consists  of  two  distinct  modules  where  airfield  aggregation 
and  route  selection  are  performed.  The  first  module  aggregates  origin  and  destination 
airfields  into  selected  centroid  airfields.  The  second  module  generates  the  route  set  which 
will  be  used  in  the  optimization  model. 

A.  CENTROID  MODEL 

As  outlined  in  the  introduction,  aggregation  of  airfields  may  lead  to  infeasibilities  in 
the  solution  produced  by  the  resulting  strategic  airlift  model  with  respect  to  borrowing 
airfield  capacity  and  with  respect  to  an  aircraft’s  ability  to  fly  certain  flight  legs  due  to 
altered  inter-airfield  distances.  It  is  primarily  for  the  latter  reason  that  we  propose  a 
centroid  aggregation  model  that  has  a  primary  objective  of  minimizing  the  maximum 
centroid-satellite  distance.  A  pure  minimax  model  will  typically  have  many  multiple 
optimal  solutions,  but  some  optima  are  usually  better  than  others  with  respect  to  criteria 
ignored  by  the  minimax  objective  (e  g.,  the  satellite-centroid  for  second  furthest  airfield). 

As  a  result  we  incorporate  a  secondary  minisum  objective  term  in  the  centroid  selection 
model  to  resolve  multiple  optimal  solutions.  Aggregation  is  not  performed  for  enroute 
airfields,  but  is  performed  for  origins  and  destinations  separately. 

The  preprocessing  step  is  based  on  a  set  of  raw  data  provided  by  various  sources 
such  as  the  TPFDD  and  airfield  infrastructure  data.  Due  to  some  policy  constraints  or  data 
inaccuricies,  there  are  occasionally  some  triangle  inequality  violations  for  inter-airfield 
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distances  As  a  data  integrity  check  we  ensure  that  given  any  two  points  X  ~ 

and  Y  =  (>’,  the  great  circle  distance  between  X and  7,  which  we  denote  by  \X -  V\  has 

the  following  properties  for  all  X  and  V: 

1  Nonnegativity  \X  -  K|  >  0 

2,  Symmetry  |A"  -  K|  =  [7  -  A"| 

3.  Triangle  Inequality  \X -Y\<\X  -  Z\  +  \Z -Y\  for  any  Z 

For  the  given  data  set,  if  the  triangle  inequality  violation  occurs, 

\X -Y\>\X -Z\-^\Z-Y\  foranyA:r,  Z 
an  algorithm  is  used  to  replace  this  inequality  with  the  following  equality. 

|A-r|  =  |A'-z|  +  |z-r| 


1.  Mathematical  Formulation 

a.  Index  Sets 

i  el  :  Original  airfields,  e  g.,  Travis  AFB 

y  e./  :  Candidate  airfields  that  may  be  centroids 

h.  Data 

d,j :  Great  circle  distances  between  airfield  /  and  centroid  j 
k  :  Number  of  origin  (or  destination)  airfields  that  will  be  used  in 
the  Throughput  II  model 
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Decision  Variables 


if  airfield  is  assigned  to  centroid  j 


X  - 

Q  otherwise 


1  if  candidate  airfield  j  is  selected  as  a  centroid 
0  otherwise 


Minimax  Model 


min  max 

s.t.  x_^<y^ 


V{i,j)sIxJ 


V/  e/ 


(2.1) 


g{0,1}  V{ij)eIxJ 
yj  ^  ^ 

The  integer  program  in  (2.1)  is  a  variation  of  Schrage  (1991)  that  may  be 
converted  to  a  linear  integer  programming  problem  by  introducing  a  new  decision  variable 


s.t.  Z>d,jX,j 

Z->'y  =  ^ 

JiJ 


V(/,7)  g  /  X  J 
W{ij)GlxJ 


y  i  si 


y{ij)sixj 
yj  sj 


(2.2) 
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This  model  minimizes  the  maximum  distance  between  centroids  and  their 


assigned  airfields  while  ensuring  that  all  airfields  are  assigned  to  only  one  centroid  and  the 
total  number  of  centroids  equals  a  prespecified  value.  For  reasons  detailed  above,  we  use  a 
secondary  minisum  objective  to  select  among  multiple  optimal  solutions.  This  can  be 
stated  mathematically. 


li.l  I(-.J 

'  ^ 

Al 

V(/,y)G/x  J 

X  <  y 

U  ■'I 

Z>’./ 

nJ 

V(/,_/)  g/  X  ./ 

II 

y  i  G  1 

x,  e(0,1) 

V(/, y)  gI  X  ,/ 

|o.'} 

Vy'  gJ 

where  f  >  0  is  a  sufficiently  small  prespecified  constant.  Model  (2.3)  contains  a  large 
number  of  binary  variables  that  can  quickly  result  in  unreasonably  long  solution  times, 
even  for  modest  values  of  |/|  and  \  j\ .  As  a  result,  we  utilize  a  solution  technique  that 
solves  a  sequence  of  a  set  covering  location  models  to  find  the  optimal  solution  to  (2.3). 

e.  Set  Covering  Model 

The  covering  problem  involves  selecting  the  minimum  number  of  centroids 
so  that  each  of  the  original  airfields  is  within  a  specified  distance  of  some  centroid  (i.e., 
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covered  by  that  centroid)  (Francis  et.  al.,  1992).  The  set  covering  model  can  be  stated  by 
introducing  the  following  addition  notation, 

(l)Data 

r  Maximum  allowed  distance  between  a  centroid  tind  its  assigned 
airfields 


(2)  Index  Set 


n,  =  {  <  r 


} 


set  of  centroids  that  may  used  to  cover  airfield  I 


J(-J 

^  Y.yj  V  /  e/ 


(2.4) 


The  key  connection  between  the  set  covering  model  (2.4)  and  minimax -minisum  model 
(2.3)  is  the  following.  If  Zj,  ,  the  optimal  objective  value  to  (2.4)  for  a  specific  value  of  r, 

exceeds  k  ( the  centroid  budget  in  (2.3))  then  we  know  r*  3,  >  r.  Similarly  if  z*^^^  <  k 

then  we  know  ^  Our  goal,  then,  is  to  find  the  smallest  value  of  r  for  which  (2.4) 

yields  an  objective  of  at  most  k.  Such  an  approach  proves  to  be  computationally  efficient 
because: 

(i)  the  number  of  binary  decision  variables  in  (2.4)  is  |,/[  while  the  nun:d)er  of  binary 
variables  in  (2.3)  is  |/|  •  |./|  +  |y|  and 
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(ii)  empirical  evidence  has  shown  that  a  linear  relaxation  of  (2.4)  often  provides  a  “near¬ 
integer”  solution. 

The  algorithm  for  solving  (2.3)  via  a  sequence  of  set  covering  models  as  defined  in 
(2.4)  is  shown  in  Figure  2.1.  The  main  part  of  the  algorithm  solves  (2.3)  vidth  £■  =  0,  and  a 
postprocessing  step  solves  an  appropriately  defined  minisum  model  to  resolve  multiple 
optimal  minimax  solutions  according  to  the  secondary  minisum  objective. 


Figure  2. 1  Set  Covering  Location  Algorithm 


2. 


The  post-processing  minisum  model  may  be  stated  : 


Index  Set 


n  =  {(/,., 

VI 

min 

a  7  ji  J 

V! 

v(/,  ./)  gO 

II 

II 

< 

m 

y,  £!».!) 

y  j  eJ 

<1 

y{‘j) 

Note  that  due  to  the 

linearity  of  the  objective  the  binary  constraints  on  Xy  can  be 

relaxed  in  (2.5)  (Schrage,  1991).  This  mixed  integer  program  selects  a  solution  from  the 
set  of  optimal  solutions  to  (2  .1),  because  centroid  assignments  are  approximately 
restricted  by  the  set  Q .  General  Algebraic  Modeling  System  (GAMS)  (Brooke  et.  al., 
1992)  code  for  solving  model  (2.3)  via  the  algorithm  of  figure  2. 1  is  in  Appendix  B. 
Computational  results  for  this  and  alternative  approaches  are  discussed  in  Chapter  V 

B.  ROUTE  SELECTION  MODEL 

The  primary  goal  of  strategic  airlift  planning  is  to  satisfy  Commander-in-Chief 
Central  Command  (CrNCCENT)’s  requirements  by  employing  airlift  resources  effectively. 
To  meet  the  goal,  planners  have  to  consider  the  entire  airlift  system  and  its  interrelated 
parts.  Each  enroute  stop  and  destination  airfield  must  have  adequate  runways,  taxiways. 
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ramps  and  support  facilities  Where  crew  changes  are  required,  a  sufficient  number  of 
qualified  and  properly  rested  aircrews  must  be  available.  Normal  crew  duty  time 
limitations  (unless  otherwise  specified  in  the  governing  OPORD/OPLAN)  are  16  hours  for 
a  basic  crew  and  24  hours  for  an  augmented  crew. 

The  route  selection  model  may  be  viewed  as  a  search  process  or  more  specifically 
as  a  traversal  of  a  tree  structure  in  which  each  node  represents  an  airfield  and  each  arc 
represents  a  flight  leg  between  the  airfields  represented  by  the  nodes.  The  search  process 
must  find  paths  through  the  tree  that  connect  the  origin  with  the  destination  and  obey 
given  physical  constraints  for  each  aircraft  type  such  as  the  appropriate  airfield  types  it 
may  utilize,  the  aircraft  minimum  payload  ferry  range,  crew  stages  etc.  Because  of  the 
extremely  large  number  of  possible  routes  that  could  be  taken,  we  need  to  bound  the 
number  of  enroute  stops.  Governing  our  search  is  a  set  of  rules  that  defines  the  legal 
routes  we  can  follow.  These  rules  include  the  crew  stages,  enroute  airfields,  runway 
limitations  in  airfields  and  the  type  of  airfield  (e  g.,  military,  civilian,  joint).  Suppose  that 
we  want  to  find  some  routes  from  node  A  to  node  G  through  a  network  of  airfields,  such 
as  shown  in  Figure  2.2. 


Figure  2  2  A  Basic  Route  Selection  Problem 
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Finding  a  route  involves  six  kinds  of  efforts: 

•  Each  leg  of  the  flight  should  be  less  than  the  minimum  payload  ferry  range  of  the 
corresponding  aircraft, 

•  The  cumulative  time  including  ground  time  spent  at  each  stop,  without  a  crew  change 
should  be  less  than  24  hours, 

•  Every  airfield  on  the  route  should  be  compatible  with  the  aircraft, 

•  The  length  of  the  route  should  be  less  than  the  upper  bound^  which  is  defined  by  the 
user, 

•  All  airfields  can  be  used  by  all  types  of  aircraft  as  an  APOD  or  an  APOE  but  only 
specified  aircraft  as  an  enroute, 

•  The  distance  between  two  adjacent  enroute  airfields  should  be  greater  than  a 
prespecified  distance  defined  by  the  user. 

The  most  obvious  way  to  find  a  solution  is  to  devise  a  bookkeeping  scheme  that 
allows  an  orderly  representation  of  all  possible  paths.  It  is  useful  to  note  that  the 
bookkeeping  scheme  must  not  allow  a  node  to  be  revisited  (i.e.,  a  cycle)  in  the  network. 
With  cyclic  paths  eliminated,  networks  are  equivalent  to  trees.  The  tree  shown  in  Figure 
2  3  is  made  from  the  network  in  Figure  2.2  by  using  breadth  first  search  and  eliminates 
certain  routes  (e  g.,  ADBG)  by  using  physical  and  strategic  constraints.  Pseudo  code  for 
the  route  generation  algoritm  is  described  below  and  additional  details  are  contained  in 
Appendix  D 

■’  Upper  bound  is  defined  by  a  parameter  times  the  great  circle  distance  between  origin  and  destination. 
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1.  Pseudo  Code  of  Route  Generation  Model 

read  input  data  files 

compute  the  great  circle  distances  between  airfields  and  add  a  random  number  between  0-1^ 
for  dclivcr>'/rccover>^  routes 
for  each  aircraft 

determine  the  adjacency  relation  between  airfields 

if  (distance  <  ferry  range  )  and  (airfield  i  <>  airfield2  )  and 
(aircraft  is  compatible  with  airfield  )  then 
adjacent  =1.0 
otherwise 

adjacent  =  0.0 

end  if 

route  generation  algorithm 

put  the  origin  as  first  stage  in  the  tree 

LOOP  :  stage  2  to  stage  5 

LOOP  :  if  eof  (previous  stage)  then  exit  (of  LOOP) 
node:=  first(previous  stage) 
if  node  destination  then  begin 

select  nodes  which  are  accessible  from  node  due  to 
given  constraints 

exclude  any  node  which  is  found  in  the  past  route  list 
add  rest  of  the  nodes  to  next  stage 


^  If  the  geographic  location  values  are  not  available  based  on  security  reason,  program  can  just  read  the 
distance  matrix  provided  by  user.  By  adding  a  small  random  number,  the  program  perturbs  the  data  in 
order  to  facilitate  breaking  ties  and  having  routes  uniquely  identified  by  their  length. 
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end 

goto  LOOP 
increment  stage  level 
end  LOOP. 

LOOP  start  from  the  last  stage  to  first  stage 
find  the  destination  node 
follow  past  route  list  through  origin 
record  the  route 
end  LOOP, 
feasibility  checks 

if  route  length  >  upper  bound  then  delete  the  route 
if  route  is  infeasible  due  to  crew  stages  then  delete  the  route 
manual  screening  of  the  routes  by  user 

enumerate  routes  based  on  their  length  and  store  for  each  aircraft 

end  for 

manual  route  adding  by  user 

end  for 

output  the  route  set. 
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IV. 


POSTPROCESSING  MODULES 


The  final  output  of  the  combined  centroid  and  route  selection  preprocessing  step  is 
a  set  of  candidate  routes  that  can  be  used  as  input  data  for  Throughput  II.  However  the 
airfield  aggregation  and  subsequent  route  selection  process  may  still  generate  a  large 
number  of  routes,  and  it  may  not  be  computationally  practical  to  include  all  the  candidate 
routes  in  the  Throughput  II  model.  Thus,  we  utilize  a  method  to  convert  the  large 
problem  into  smaller  problems  of  manageable  size.  The  approach  to  this  reduction  is  to 
solve  a  linear  program  optimized  over  a  subset  of  its  columns,  the  result  of  which  may  be 
suboptimal  Based  on  ideas  from  LP  sensitivity  analysis  we  can  test  whether  the  set  of 
decision  variables  associated  with  a  certain  route,  not  in  the  original  subset,  is  attractive.  If 
it  is  attractive,  it  can  be  put  into  the  basis  as  a  new  column.  Moreover,  we  can  examine 
why  this  route  might  be  attractive  and,  if  necessary,  modify  the  rules  that  select  the  initial 
set  of  routes  in  the  heuristic  route  selection  model.  The  Throughput  II  model  may  be 
regarded  as  an  example  of  the  following  LP  model  : 

min  cx  (3  1) 

St  Ax  =  h 

X  >  0 

and  its  dual 

mdiX  nh  (3.2) 

s.l  n  A  <  c 

n  unrestricted . 
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The  optimality  conditions  for  a  linear  program  state  that  necessary  and  sufficient 
conditions  for  x*  and  n'  to  be  an  optimal  primal-dual  pair  are, 


1 

Primal  Feasibility 

Ax  h 

2. 

Dual  Feasibility 

71  A  <c 

Strong  Duality 

cx  =  71  h 

In  standard  column  generation  techniques,  new  columns  are  defined  with  respect 
to  the  original  constraint  set.  Throughput  II,  however,  does  not  fall  into  this  framework. 
Instead,  associated  with  each  new  route-aircraft  combination  at  any  point  in  time  is  a  set 
of  decision  variables  that  appear  in  the  original  constraints  plus  an  additional  set  of 
constraints  that  contain  only  the  new  decision  variables.  For  example,  one  of  the  new 
constraints  specifies  that  the  amount  of  cargo  that  can  flow  along  a  new  route  is  restricted 
by  the  capacity  of  aircraft  flying  this  route.  We  denote  the  set  of  variables  associated  with 
the  new  routes >’  and  express  their  inclusion  in  model  (3.1)  via, 

min  cx  +  ^  (3.3) 

sA.  Ax  A  h'y  -  h 
Dy^O 
x.y>0 

The  fact  that  the  right-hand-side  of  the  new  set  of  constraints  is  zero  reflects  our 
assumption  that  the  new  LP  (3.3)  is  feasible  \^y  0  ,  i.e.,  we  are  not  required  to  fly  the 
new  routes.  The  dual  of  (3.3)  is. 
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max  nh 
s.t  n  A  <  c 


(3  4) 


n  F  y  D  <  f 

As  a  result  of  our  non-standard  form  the  usual  optimality  test  regarding  the  sign  of 
the  reduced  cost  of  a  new  column  no  longer  directly  applies.  However,  the  following 
theorem  provides  an  analogous  optimality  test  for  the  framework  outlined  above. 

Theorem  3.1 

Let  X*  and  n'  be  an  optimal  primal-dual  pair  for  (3.1)  and  (3.2).  If  nFj  -  fj  <0  \/j , 
then  (x* ,  0)  and  {n\Qi)  form  an  optimal  primal-dual  pair  for  (3.3)  and  (3.4). 

Proof  of  Theorem  3.1 : 

We  will  verify  the  Karush  -  Kuhn  -  Tucker  (KKT)  Optimality  Conditions  of 

1  Primal  Feasibility 

2  Strong  Duality 

3  Dual  Feasibility 

of  the  primal-dual  solution  pair  (x",y'’)  =  (x*,0)  and  (7r",y")  =  (tt*  ,0)  in 
(3.3)  and  (3.4). 
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1.  Primal  Feasibility  : 

Ax"  +  Fy"  =  h 
Dy"  =  0 
x'\y"  >  0 

Now, 

Ax'  +FQ-h  since  Ax*  b  and 
DO  =  0 

(x*  ,0)  >  0  also  follows  from  feasibility  of  x*  in  (3.1). 

2.  Strong  Duality  : 

cx"  +  fy"  -  n" h  +  y  "0 
cx*  +  /  0  =  K*  h  +  yO 

+  *  L 

CX  -no 

follows  from  strong  duality  of  x*  and  n*  in  (3. 1)  and  (3.2). 

3.  Dual  Feasibility  : 

Tc"  A  ^  -  Cj  =  TT*  A  j  -  Cj  <0  follows  from  dual  feasibilty  in  (3.2). 
n"F^-y"D^-f^^n*F,-f^<0 
follows  from  the  optimality  test  hypothesis.  QED 

From  Theorem  3  . 1  we  know  that  if  we  solve  Throughput  II  and  subsequently 
determine  that  the  reduced  costs  of  variables  associated  with  routes  not  explicitly  included 
in  the  original  model  pass  the  optimality  test  k*  F ^  -  /^  <  0  then  the  current  solution  is 
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optimal  with  respect  to  all  routes  If  the  optimality  test  is  not  satisfied  then  we  can  append 
the  variables  associated  with  attractive  routes  and  re-solve  the  corresponding  Throughput 
11  model  To  implement  the  column  generation  technique  in  Throughput  II,  note  that  the 
variables  associated  with  the  new  routes  are  X,  Y,  TONSUE,  and  TPAX.  Thus,  summing 
the  dual  variables  associated  with  constraints  that  these  variables  appear  in  computes  the 
“;r/*  ”  terms  Additionaly,  only  the  TONSUE  and  TPAX  variables  have  objective  function 
coefficients  (for  late  penalty  computation).  Consequently,  only  these  columns  have 
terms  associated  The  complete  GAMS  code  is  included  in  Appendix  C. 

It  is  desirable  to  prevent,  or  at  least  recognize,  when  Throughput  II  generates 
solutions  that  are  infeasible  with  respect  to  the  original,  disaggregated  network.  One 
example  of  such  a  violation  may  occur  when  a  unit’s  origin  airfield  (and  hence  the  unit)  is 
aggregated  with  a  centroid  closer  to  the  destination  airfield;  in  this  case,  the  origin- 
destination  distance  may  be  greater  than  the  selected  aircraft’s  minimum  payload  range, 
even  though  the  centroid-destination  distance  was  not  longer  than  the  minimum  payload 
ferry  range  In  reality,  therefore,  this  aircraft  wold  have  to  use  a  different  route  or  carry  a 
lighter  load.  Another  potential  infeasibility  can  arise  when  Throughput  II  inappropriately 
“borrow”  airfield  capacity  between  airfields  that  have  been  aggregated  into  a  single 
airfield.  A  promising  and  important  area  for  future  research  is  to  implement  a 
disaggregation  procedure  that  maps  aircraft  missions  to  the  original  network.  With  such  a 
procedure  in  place,  it  will  be  possible  to  recognize  and  take  steps  to  avoid  the 
infeasibilities  described  above. 
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V. 


RESULTS  AND  SUMMARY 


This  chapter  reports  the  results  for  model  instances  based  on  two  different  data 
sets  provided  by  AFSAA  in  the  summer  and  fall  of  1995.  TPFDD  data  must  be 
aggregated  for  use  in  Throughput  II.  This  aggregation  procedure  is  implemented  in  the 
Pascal  programming  language.  It  takes  the  input  TPFDD,  and  then  performs  aggregation 
based  on  a  strategy  determined  by  several  user-specified  parameters.  The  flow  diagram  of 
the  aggregation  algorithm  is  in  Appendix  D  and  may  be  summaiized  as  follows: 

Step  1 :  Aggregate  units  that  are  to  be  shipped  between  the  same  origin-destination 
pair  and  have  identical  delivery  time  windows  based  on  the  an^lable-load-date  (ALD) 
and  required-delivery-date  (RDD) . 

Step  2:  Eliminate  units  that  have  very  small  movement  requirements  (e.g.,  one 
passenger). 

Step  3 :  Aggregate  units  that  are  to  be  shipped  between  the  same  origin-destination 
pair  and  have  similar  RDD  values.  The  definition  of  “similar”  is  allowed  to  depend  on 
when  over  the  deployment  period  the  units  should  be  delivered.  Greater  latitude  is 
permitted  near  the  end  of  the  scenario  compared  to  the  early  days  of  the  deployment.  This 
is  implemented  by  the  user  specifying  parameters  Ri,  R2,  R3,  partitioning  the  planning 
period  into  segments  with  corresponding  tolerance  values  r/,  rj,  where, 

0<r,<r,<r, 
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If  two  units  with  the  same  origin  and  destination  have  RDD’s  in  the  same  segment 
and  if  the  difference  between  the  RDD’s  is  less  than  the  corresponding  r  value,  the  two 
units  are  aggregated  This  is  depicted  in  Figure  5.1. 


\RDD}-RDD2  :  <r, 

1 

RDD!-RDD2\  <  r, 

1 

:  RDD  I  -  RDD2  \  <  r, 

1 

0 

Ri 

R2 

R3 

Figure  5  1.  Unit  Aggregation  Over  the  Time  Horizon 


After  the  aggregation  the  RDD  and  ALD  values  are  defined  as: 

RDD„,„.=  Max  (RDD1.RDD2) . 

ALD„c„  =  Weighted  Based  on  Movement  Requirements  (ALD1,ALD2) 

Step  4:  The  results  of  the  aggregation  procedure  are  then  used  in  the  centroid 
model,  followed  by  the  route  selection  procedure  which  in  turn  provides  input  data  for  the 
Throughput  1 1  optimization  model, 

A.  RESULTS 

1.  Centroid  Model 

The  data  set  that  has  been  used  for  the  centroid  models  consists  of  1 1 7  airfields  in 
the  Continental  United  States  (CONUS)  All  the  models  presented  in  Chapter  III  can  be 
used  to  select  the  centroid  airfields.  In  order  to  investigate  the  relative  merits  of  these 
models  three  criteria  have  been  defined.  The  first  one  is  the  computational  time  required  to 
solve  the  model  on  an  IBM  RS6000/590  system  using  GAMS  with  the  OSL  solver.  The 
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second  criterion  is  the  average  distance  between  the  centroid  airfields  and  their  satellite 
airfields.  This  is  an  important  issue  in  order  to  have  an  accurate  delivery  cycle  time  in  the 
Throughput  II  optimization  model.  The  last  criterion  is  the  primary  objective  value  of  the 
proposed  model  which  is  the  maximum  distance  between  the  centroid  airfields  and  their 
satellite  airfields.  The  model  that  we  develop  allows  the  user  to  predefine  certain  centroid- 
satellite  groups  and  effectively  remove  them  from  consideration  in  the  centroid  model. 

This  is  desirable  feature  because  it  allows  the  user  to  prevent  the  possibility  of  capacity 
borrowing  at  known  system  bottlenecks.  The  results  are  illustrated  in  Table  5.1. 


MODEL 

TIME 

(RS6000/590 

GAMS/OSL) 

AVERAGE 

DISTANCE 

(nm) 

MAXIMUM 

DISTANCE 

(nm) 

YOST’S  MODEL 

14  SECONDS 

326 

794 

MINISUM  MODEL 

10  MINUTES 

277 

717 

MINIMAX 

MODEL^ 

32  SECONDS 

370 

600 

SUGGESTED 

MODEL 

87  SECONDS 

317 

600 

Table  5.1.  Computational  Results  for  Various  Centroid  Models. 


Note  that  Yost’s  model  is  a  minisum  model  with  a  limited  number  of  airfields  that 
may  be  used  as  centroids.  As  a  result,  in  runs  quickly  but  performs  worse  than  the  other 
models  relative  to  the  other  criterion.  The  suggested  model  provides  an  attractive  balance 
between  the  minisum  model  and  the  minimax  model  with  respect  to  the  average  and 
maximum  distance  criterion  and  requires  modest  computational  effort 

This  is  the  model  that  solved  via  set  eovering  algorithm. 
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Recall  that  the  centroid  models  require  the  user  to  define  the  (maximum)  number 


m 


2. 


Route  Selection  Model 


After  we  select  the  centroid  airfields  for  origin  and  destination  airfields  via  the 
centroid  model,  the  route  generation  model  constructs  routes  based  on  user-defined  input 
parameters  It  takes  less  than  an  hour  to  set  up  the  input  files  for  the  first  run.  After  this 
preparation,  each  run  takes  approximately  five  minutes. 

The  first  input  parameter  is  denoted  Z  and  is  used  to  define  the  maximum  route 
length  between  origin-destination  pairs.  This  lambda  value  actually  consists  of  the  product 
of  three  different  parameters.  We  assign  a  lambda  value  for  each  airfield  and  aircraft.  To 
illustrate  this  process,  if  we  call  the  resultant  lambda  value  as  the  composite  lambda  value, 
then. 


A 


composite 


origin 


*  X 


destination 


*X 


aircraft 


If  we  define  each  of  these  parameters  as  1.0,  then  only  the  direct  route  (i.e.,  no 
enroute  stops)  will  be  generated  Using  these  three  different  parameters  provides  flexibility 
for  selecting  routes  based  on  the  type  of  aircraft  or  accessibility  of  different  regions  of  the 
world  with  origin/destination  dependent  lambda  values. 

The  second  parameter  of  the  route  selection  model  is  the  restriction  on  the  flight 
time  of  an  aircrew  In  this  research  we  used  24  hours  for  the  maximum  flight  time  without 
a  crew-change.  The  final  input  required  from  the  user  is  the  set  of  crew-stage  airfields  and 
enroute  airfields  that  can  be  used  during  the  airlift  operation.  Table  5.2  shows  the  number 
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of  infeasible  routes  for  various  X  values  for  the  C747  which  has  a  minimum  payload 
ferry  range  of  6500  nm  for  delivery  routes  and  9200  nm  for  recovery  routes. 


Total  Routes 

Composite 

Lambda  Value 

Total 

Feasible  Routes 

Total  Infeasible  Routes  Due 

To  Flight  Time  Restrictions 

1.041 

116 

0 

7266 

1.117 

163 

6 

Delivery  Routes 

1.161 

217 

22 

1.184 

254 

27 

1.240 

302 

43 

1.000 

32 

0 

4152 

1.139 

149 

26 

Recovery  Routes 

1.162 

233 

53 

1.185 

274 

73 

Table  5.2  Route  Selection  Model  Results 


We  have  performed  experiments  in  order  to  show  the  “quality”  of  the  routes  that 
have  been  generated  by  the  route  selection  procedure  versus  the  routes  with  the  same  data 
set  in  Lim  (1994)  The  computational  results  based  on  comparing  manually  generated 
routes  and  the  model  are  shown  in  Table  5.3, 
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Model 

Number  of 
Routes 

Objective  Function 
Value 

Route  Selection  Model 

71  routes 

13.208 

Manually  Generated 

Routes 

62  routes 

15.021 

Table  5  3  Route  Selection  Model’s  Results. 


The  route  selection  model  improves  the  routing  options  in  Throughput  II  so  that 
there  is  a  12%  reduction  in  the  weighted  sum  of  late  and  undelivered  cargo. 

3.  Postprocessing  Module 

One  important  issue  in  the  postprocessing  module  is  whether  including  additional 
routes  in  Throughput  II  will  result  in  an  improved  solution.  There  are  many  decision 
variables  associated  with  a  single  physical  route  through  the  network: 

The  route  may  be  flown  by  several  types  of  planes  and  at  different  points  in  time  in 
the  planning  period  There  are  typically  a  large  number  of  decision  variables  associated 
with  routes  not  explicitly  included  in  Throughput  II  that  have  attractive  reduced  costs  (see 
Chapter  II).  However,  many  of  these  provide  small  if  any  improvement  in  the  solution,  if 
included.  For  this  reason,  we  only  augment  the  set  of  routing  options  with  decision 
variables  whose  reduced  costs  are  significantly  attractive,  relative  to  some  tolerance.  In 
the  proposed  method.  Throughput  II  is  solved  with  respect  to  the  original  routes  and  then 
the  set  of  routes  is  augmented  and  the  model  is  re-optimized. 
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When  this  is  done,  it  is  desirable  to  initialize  the  re-optimization  with  the  basis 
from  optimal  solution  of  the  original  model.  Table  5.4  illustrates  the  savings  in 
computational  effort  that  this  technique  provides  relative  to  solving  the  larger  model  from 
scratch.  In  addition.  Table  5.4  shows  that  the  postprocessing  scheme  found  routing 
options  that  significantly  decreased  the  weighted  sum  of  late  and  undelivered  cargo. 


Number  of 
Routes 

Model 

Objective 
Function  Value 

Running  Time 
(RS6000) 

71  routes 

Throughput  11 

13.208 

1 0  minutes 

117  routes 

Throughput  II 

8.0292 

60  minutes 

117  routes 

Proposed  Model 

8.0292 

38.75  minutes 

Table  5.4  Computational  Results  of  Sensitivity  Analysis 

B.  RECOMMENDATIONS  AND  FUTURE  WORK 

The  techniques  proposed  in  this  thesis  have  been  recently  applied  to  a  large 
TPFDD  provided  by  AFSAA  in  August  1995  for  a  C17/C747  fleet  mix  study.  Due  to  the 
enormous  size  of  the  data,  which  has  approximately  14500  movement  requirements,  the 
preprocessing  modules  for  centroid  selection  and  route  reduction  were  an  invaluable  tool 
for  quickly  generating  an  optimization  model  of  tractable  size.  This  modeling  and  data 
reduction  exercise  also  revealed  that  certain  additional  features  could  further  streamline 
the  process  for  future  analyses: 

•  User  interface.  Current  aggregation  and  route  selection  modules  are  implemented  as  a 
Disk  Operating  System  (DOS)  application.  Because  of  the  structure  of  this  operating 
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system,  the  ease  of  use  for  an  analyst  would  be  greatly  enhanced  by  implementing  this 
software  as  a  Windows  application  with  user-friendly  interfaces. 

•  Data  Interface:  The  current  model  uses  text  files  that  must  be  manually  entered  by  the 
analyst  Because  this  is  a  time  consuming  procedure,  a  data  interface  that  permits  easy 
modifications  while  ensuring  data  integrity  would  significantly  decrease  the  model 
generation  time. 

•  Route  Generation  Module:  As  indicated  above,  there  is  a  set  of  decision  variables 
associated  with  each  route.  It  would  be  useful  to  develop  a  theory  that  utilizes  reduced 
cost  information  from  individual  decision  variables  to  indicate  when  the  inclusion  the 
set  of  decision  variables  will  improve  (or  be  likely  to  improve)  the  current  solution. 

•  Feasibility  Test  Postprocessing  Module:  As  discussed  at  the  end  of  Chapter  IV, 
deriving  a  procedure  that  maps  Throughput  II  solutions  to  the  original  disaggregated 
network  would  prove  valuable  with  respect  to  testing  the  feasibility  of  the  proposed 
solution 

As  shown  above,  much  work  remains  to  realize  the  full  potential  of  pre  and  post¬ 
processing  of  the  Throughput  II  model.  The  contributions  of  this  thesis  form  a  solid 
foundation  for  these  improvments. 
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APPENDIX  A.  THROUGHPUT  II  MODEL 


The  following  is  a  brief  summary  of  Throughput  II  (Morton  et.  ai,  1995), 

A.  INDICES 

14  indexes  units,  e.g.,  82nd  Airborne 

a  indexes  aircraft  types,  e  g.,  C5,  Cl 41 

t,  t  ’  index  time  periods 

h  indexes  all  airfields  (origins,  enroutes  and  destinations) 

/  indexes  origin  airfields 

k  indexes  destination  airfields 

r  indexes  routes 

B.  INDEX  SETS 

1.  Airfield  Index  Sets 

B  set  of  available  airfields 

origin  airfields 

K^B  destination  airfields 


2.  Aircraft  Index  Sets 

A  set  of  available  aircraft  types 

A  huik^A  aircraft  capable  of  hauling  bulk  size  cargo 
A  bulk  aircraft  capable  of  hauling  over-sized  cargo 
Aout^over  aircraft  capable  of  hauling  out-sized  cargo 

Bulk  cargo  is  palletized  on  88  x  1 08  inch  platforms,  which  can  fit  on  any  plane. 

Over-sized  cargo  is  non-palletized  rolling  stock;  it  is  larger  than  bulk  cargo  and  can  fit  on 

aC141,C5orC17,  Out-sized  cargo  is  very  large  non-palletized  cargo  that  can  fit  into  a 

C5  or  C17  but  not  a  C141. 
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Route  Index  Sets 


3. 

R  set  of  available  routes 

permissible  routes  for  aircraft  a 

Rahdla  permissible  routes  for  aircraft  a  that  use  airfield  b 

Ruik^ia  permissible  routes  for  aircraft  a  that  have  origin  /  and  destination  k 

l)R,di  delivery  routes  that  originate  from  origin  / 

RRk^R  recovery  routes  that  originate  from  destination  k 

4.  Time  Index  Set 

7  set  of  time  periods 

Tuar^T  possible  launch  times  of  sorties  for  unit  u  using  aircraft  a  and  route  r 

The  set  7,,„r  covers  the  allowed  time  window  for  unit  u,  which  starts  on  the  unit’s 
available-to-load  date  and  ends  on  the  unit’s  required  delivery  date,  plus  some  extra  time 
up  to  the  maximum  allowed  lateness  for  the  unit. 

C.  GIVEN  DATA 

1.  Movement  Requirements  Data 

MoveRAXuik  Troop  movement  requirement  for  unit  u  from  origin  /  to 
destination  k 

MoveUEuik  Equipment  movement  requirement  in  short  tons  (stons)  for  unit  u 
from  origin  /  to  destination  k 

ProHulku  Proportion  of  unit  u  cargo  that  is  bulk-sized 

ProOver,,  Proportion  of  unit  u  cargo  that  is  over-sized 

ProOup  Proportion  of  unit  u  cargo  that  is  out-sized 

2.  Penalty  Data 

LatePenUKu  Lateness  penalty  (per  ston  per  day)  for  unit  u  equipment 
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LaiePenPAXu  Lateness  penalty  (per  soldier  per  day)  for  unit  u  troops 

NoGoPenUh,,  Non-delivery  penalty  (per  ston)  for  unit  u  equipment 

NoGoPenPAXu  Non-delivery  penalty  (per  soldier)  for  unit  u  troops 

MaxLate  Maximum  allowed  lateness  (in  days)  for  deliveiy 

Preserve^,  Penalty  (small  artificial  cost)  for  keeping  aircraft  a  in  mobility 
system  at  time  t 

3.  Cargo  Data 

IJESqPtu  Average  cargo  floor  space  (in  sq,  ft.)  per  ston  of  unit  u  equipment 

PAXW/u  Average  weight  of  a  unit  u  soldier  inclusive  of  personal 

equipment 

4.  Aircraft  Data 

Number  of  aircraft  of  type  a  that  become  available  at  time  / 

Maximum  troop  carriage  capacity  of  aircraft  a 

Average  cargo  space  (in  sq.  ft.)  consumed  by  a  unit  u  soldier  fro 
aircraft  a 

Cargo  floor  space  (in  sq.  ft.)  of  aircraft  a 

Cargo  space  loading  efficiency  (<1)  for  aircraft  a.  This  accounts 
for  the  fact  that  it  is  not  possible  in  practice  to  fully  utilize  the 
cargo  space. 

Established  utilization  rate  (flying  hours  per  aircraft  per  day)  for 
aircraft  a 

5.  Airfield  Data 

MOGC.aph  Aircraft  capacity  (in  narrow-body  equivalents)  at  airfield  b 


Supply^, 

MaxPAXa 

PAXSqlX, 

ACSqFta 

LoadEffu 

llRatea 
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D. 


MOGReqah  Conversion  factor  to  narrow-body  equivalents  for  one  aircraft  of 
type  a  at  airfield  h 

MOGEff  MOG  efficiency  factor  (<1),  to  account  for  the  fact  that  it  is 

impossible  to  fully  utilize  available  MOG  capacity  due  to 
randomness  of  ground  times 


6.  Aircraft  Route  Performance  Data 


MaxLoadar 

Maximum  payload  (in  stons)  for  aircraft  a  flying  route  r 

GTimeahr 

Aircraft  ground  time  (due  to  onload  or  offload  of  cargo, 
refueling,  maintenance,  etc  )  needed  for  aircraft  a  at  airfield  b  on 
route  r 

DTimeahr 

Cumulative  time  (flight  time  plus  ground  time)  taken  by  aircraft  a 
to  reach  airfield  h  along  route  r 

EltTimear 

Total  flying  hours  consumed  by  aircraft  a  on  route  r 

CTimear 

Cumulative  time  (flight  time  plus  ground  time)  taken  by  aircraft  a 
on  route  r 

DaysLate^an 

Number  of  days  late  unit  m’s  requirement  would  be  if  delivered  by 
aircraft  a  via  route  r  with  mission  start  time  t 

DECISION  VARIABLES 

1.  Sortie  Variables 

Xum  Number  of  aircraft  a  that  airlift  unit  u  via  route  r  with  mission 

start  time  during  period  t 

Y„ri  Number  of  aircraft  a  that  recover  from  a  destination  airfield  via 

route  r  with  start  time  during  period  / 
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2.  Aircraft  Allocation  and  De-allocation  Variables 

AHota,,  Number  of  aircraft  a  that  become  available  at  time  t  that  are 

allocated  to  origin  / 

Releasea,,  Number  of  aircraft  a  available  at  origin  /  in  time  i  that  are  not 
scheduled  for  any  flights  from  then  on 

3.  Aircraft  Inventory  Variables 

Number  of  aircraft  a  inventoried  at  origin  /  at  time  t 
Number  of  aircraft  a  inventoried  at  destination  k  at  time  t 

NPlaneSat  Number  of  aircraft  a  in  the  air  mobility  system  at  time  t 

4.  Airlift  Quantity  Variables 

TonsUEuart  Total  stons  of  unit  u  equipment  airlifted  by  aircraft  a  via  route  r 
with  mission  start  time  during  period  t 

Rf  AXuart  Total  number  of  unit  u  troops  airlifted  by  aircraft  a  via  route  r 

with  mission  start  time  during  period  t 

5.  Elastic  (Nondelivery)  Variables 

IJENoGouik  Total  stons  of  unit  u  equipment  with  origin  /  and  destination  k 
that  is  not  delivered  in  the  prescribed  time  frame 

PAXNoGou,k  Number  of  unit  u  troops  with  origin  /  and  destination  k  who  are 
not  delivered  in  the  prescribed  time  frame 
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OBJECTIVE 


minimize 


Z  Z  Z  Z  I^alePenUK^  *  DaysLate^^^,  *  TonsUE^^, 


u  a  rk:R„  t(zT^ 


+Z  Z  Z  Z  I^^atePenPAX^  *DaysLaie^^,  *  TPAX^^^, 

u  a  rGR^tC:T^ 

+Z  Z  Z  ( NoGoPenUE^  *  UENoGo^,,  +  NoGoPenPAX^  *  PAXNoGo^,, ) 

u  i  k 

^  Preserve *  NPlanes^^ 

u  I 

The  objective  function  minimizes  the  total  weighted  penalties  incurred  for  late 
deliveries  and  non-deliveries  The  model’s  secondary  objective  is  to  choose  a  feasible 
solution  that  maximizes  unused  aircraft. 

F.  CONSTRAINTS 


X  2;  ■^ronsUE^*UENoGo,,=MoveUE,,. 


'  ^^buik  ^^atk 


\/  u,i,k:  MoveUE^,^.  >  0 


Z  Z  'ZTomUE.„*UENoGo,^>ProO„l.*MowUE,^, 


‘-I rGR,„t  ft  7’^ 


V  MoveUE^^!^  >  0 


I  Z  z 


Ui  ri  R,„t  p  7',, 


[ProOver^  +  ProOul^)*  MovelJE^^i^ ,  V  u,i,k .  MoveUE^^^.  >  0 


Z  Z  Z  77MX  +PaxNoGo  ,  =MovePAX  ^ , 

uart  uik  uik 


V  uj,k.  MovePAX >  0 
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y  y  X  ,+//,+  Release  ,  = 

/  .  /  ^  uarl  ait  ait 

r  n  DR, 

H  ,  .  +  Allot  ,  +  y  y  K  ,,  ,  V  a,i,t 

aij  1  ait  /  1  /  j  art  ’  ’  ’ 


n  t'  CTimear  1  t 


I  Kr,  +  W'-  =  ,  +  Z  Z  V  a,*,/ 


u  rsR^ 


t'-ilCTime^  ]  / 


V  aj 


NPIanes,,  -j^J^Re/ease„„.,  V  a,t 


ZZZ'f-*A-„,  +  ZZ'^-^n,- +ZZ/^. 

ryR,,f'  I  u  riR,,t'  =  ]  i  I ’-I 


(A.  10) 


<Z^l'lattes^n 


V  aj 


where 


K, . = 


/-/'  +  !,  if  t'  <t  <t'  +  CTime^^-\ 
CTime  if  t  >t'  +  CTime^^  -  1 


TPAX^^^,<MaxPAX^*X^^^„ 


\/ u,a,r,t.t  gT  (A.11) 


TomUE^^,  +  PAXWt*  TPAX,,^,  <  MaxLoad X , 


(A.  12) 


V  u,a,rj.t  er„. 


PAXSqFt^*TPAX^^^,  +  ilKSqFt^*  TonsUE^^^,  < 

A CSqFf^  *  LoadFg^  *  X^^^, ,  M  u, a,  rj:  t 


(A.  13) 


(A.  14) 


Z  Z  Z  ^  YaYj  ar*  ^ ar,  ^ 


V  I  /Rale  *  N Planes  , ,  V  a 

/  j  a  at  ^ 


ZII 


DJime^^  ]  / 


+Z  Z  2;(W007fe^^.Gr/mc*  /24).i; 

a  nR^  t'  ^{DTime^,]-t 


<  MOGEjf*  MOGCmp,  , 


yhj 


(A.  15) 


A,  2  Demand  satisfaction  constraints  for  all  classes  of  cargo 

A. 3  Demand  satisfaction  constraints  for  out-sized  cargo 

A. 4  Demand  satisfaction  constraint  for  over-sized  cargo 

A.  5  Demand  satisfaction  for  troops 

A. 6  Aircraft  balance  constraints  at  origin  airfields 

A. 7  Aircraft  balance  constraints  at  destination  airfields 

A.  8  Aircraft  balance  constraints  for  allocations  to  origins 

A.  9  Aircraft  balance  constraints  accounting  for  allocations  and  releases 

A,  1 0  Cumulative  aircraft  balance  constraints 

A.  1 1  Troop  carriage  capacity  constraints 

A.  1 2  Maximum  payload  constraints 

A  1 3  Cargo  floor  space  constraints 

A.  14  Aircraft  utilization  constraints 

A  1 5  Airfield  MOG  constraints 
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APPENDIX  B.  GAMS  CODES  FOR  CENTROID  MODEL 


STITLE  SETCOVERING 

SSTITLE  Yasin  TURKER  29JUL1995 

♦ - GAMS  AND  DOLLAR  CONTROL  OPTIONS - 

SOFFUPPER  OFFSYMLIST  OFFSYMXREF  INLINECOM  {  } 

SONMIXED 

OPTIONS 

LIMCOL  =  0  ,  LIMROW  =  0  ,  SOLPRINT  =  OFF  ,  DECIMALS  =  2 
RESLIM  =  18()()(),  ITERLIM  =  400000,  OPTCR  =0.1,  SEED  =  3141 
MIP  =  OSL 


*includc  centroid  model  for  origins 
SINCLUDE  origset.gms 
SETS  INDEX  /I  *100/ 

AF  airfields 
SINCLUDE  airficld.scl 
1(AF)  origins 
SINCLUDE  origin  l  .txl 

ALIAS  (I.IP.J); 

AL1AS(AF,AFP,AFPP); 

PARAMETER  DIS  (AF,AFP) 

SINCLUDE  disll  .txt 

SCALAR  K2  number  of  airfields  which  is  allowed  throughput  model  /4/ 
RAD  radius  of  centroid 
RBEST  best  radius  of  centroid 
UPP  upper  value  of  centroid 
LOW  lower  value  of  centroid; 

UPP=SMAX((I,J).D1S(I,J)); 

LOW=SMIN((LJ).DIS(l,J)); 

RAD=(UPP+LOW)/2; 

RBEST=UPP; 

VARIABLE 

ZORIG  maximum  number  of  centroids 
Z  lORIG  objective  fimetion  value  of  minisum  model; 

BINARY  VARIABLE 
YORIG(J)  open  centroid  j; 

EQUATIONS 

OBJ  objective  function  of  covering  problem 
OBJ  1  objective  function  of  minisum  model 

COVER(l)  covering  assignments  for  set  covering  model 
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ASSIGN  1(1. IP)  assignment  for  each  airfield  for  an  centroid 
ASS1GN2  number  of  airfield  which  is  allowed  for  throughput  model 
ASS]GN.3(1)  make  sure  cvcr>'  airfield  either  be  in  centoid  or  alone; 

♦MINIMIZE 

OBJ..  SUM(J,YORIG(J))=E=ZORlG; 

COVER(I)..  SUM(J,Y0RIG(J)$(D1S(I,J)  LE  RAD))=G=1; 

MODEL  ORIGCOV  /OBJ,COVER/; 

LOOP(INDEX$(NOT(((UPP  -  LOW)  LE  5.0)  AND  (ZORIG.L  EQ  K2))). 
RAD=(UPP+LOW)/2; 

SOLVE  ORIGCOV  USING  MIP  MINIMIZING  ZORIG; 

DISPLAY  RAD.UPP.LOW; 

UPP=RADS(ZORIG.L  LE  K2+0. 1)+UPP$(Z0R1G.L  GE  K2+0.9); 
LOW=RAD$(ZORIG.L  GE  K2+().9)+LOW$(ZORIG.L  LE  K2+0.1); 


); 

DISPLAY  YORIG.L,RAD,UPP,LOW,ZORlG.L; 

POSITIVE  VARIABLE 

XORIG(LIP)  airfield  i  belongs  to  centroid  j  ; 

X0R1G.UP(IJ)$(D1S(LJ)  LE  RAD)  =1.0; 

♦minisum  model 

OBJL.  SUM((LJ))$(DIS(LJ)  LE  RAD),DIS(LJ)*X0R1G(1,J))=E=Z10RIG; 
ASSIGN1(I,J)$(D1S(I,J)  LE  RAD)..  XOR]G(I,J)=L=YORIG(J); 

ASSIGN2..  SUM(J,YOR]G(J))=E=K2; 

ASSIGN3(I)..  SUM(J$(DIS(LJ)  LE  RAX»),XORIG(LJ))  =E=1; 

MODEL  CENTROID/OBJ  1  .ASSIGN  1 , ASSIGN2, ASSIGN3/; 

SOLVE  CENTROID  USING  MIP  MINIMIZING  ZIORIG; 

DISPLAY  YORIG.L, ZORIG  L.ZIORIG  L.XORIG.L; 


♦include  centroid  model  for  dcstinatinations 
SINCLUDE  dcstsct.gms 
SETS  K(AF)  airports 
SINCLUDE  dcstinal  txt 

ALIAS  (K,KP,L); 

UPP=SMAX((K,L),DIS(KX)); 

L0W=SM1N((K,L),DIS(KX)); 

RAD=(UPP+LOW)/2; 

SCALAR  K1  number  of  centroid  destination  airfields  is  allowed  in  Throughput  /8/; 
VARIABLE 

ZDEST  maximum  number  of  centroids 
ZIDEST  objective  function  value  of  minisum  model; 

BINARY  VARIABLE 
YDEST(L)  open  centroid  L; 
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EQUATIONS 

OBJ4  objective  function  of  covering  problem 
OBJ6  objective  function  of  minisum  model 
COVER  1(K)  covering  assignments  for  set  covering  model 
ASSIGN  1  l(K.KP)  assignment  for  each  airfield  for  an  centroid 
ASSIGN2 1  number  of  airfield  which  is  allowed  for  throughput  model 
ASSIGN!  I(K)  make  sure  every  airfield  either  be  in  centoid  or  alone; 

♦MINIMIZE 

OBJ4.,  SUM(L,YDEST(L))=E=ZDEST; 

COVERl(K)..  SUM(L,YDEST(L)$(D1S(K,L)  LE  RAD))=G=1; 

MODEL  DESTCOV/OBJ4.COVER1/; 

LOOP(lNDEX$(NOT(((UPP  -  LOW)  LE  5.0)  AND  (ZDEST.L  EQ  Kl))), 
RAD=(UPP+LOW)/2; 

SOLVE  DESTCOV  USING  MIP  MINIMIZING  ZDEST; 
UPP=RAD$(ZDEST.L  LE  Kl+0. 1)+UPP$(ZDEST.L  GT  Kl+0.9); 
LOW=RAD$(ZDEST  L  GT  Kl+0  9)+LOW$(ZDEST.L  LE  Kl+0.1); 


); 

DISPLAY  YDEST  L,RAD,UPP,LOW,ZDEST.L; 

POSITIVE  VARIABLE 

XDEST(K,KP)  airfield  i  belongs  to  centroid  L  ; 

XDEST.UP(K,L)$(DIS(K,L)  LE  RAD)  =1.0; 

♦minisum  model 

OBJ6..  SUM((K,L)  $(DIS(K,L)  LE  RAD),DIS(K,L)^XDEST(K,L))=E=ZIDEST; 
ASSIGN  1 1(K,L)$(D1S(K,L)  LE  RAD)..  XDEST(K,L)=L=YDEST(L); 
ASSIGN2L.  SUM(L.YDEST(L))=E=K1; 

ASSIGN!  1(K).  SUM(L$(DIS(K,L)  LE  RAD),XDEST(K,L))  =E=1; 

MODEL  CENTROID1/OBJ6ASS1GN1 1,ASSIGN21,ASSIGN31/; 

SOLVE  CENTROID  1  USING  MIP  MINIMIZING  ZIDEST; 

D1  SPLAY  YDEST.L,ZDEST.L,Z  1  DEST.L,XDEST.L; 


♦include  set  of  units 
SET  UNIT  units  in  TPFDD 
/ 

SINCLUDE  unitl  .txt 

/; 

SETS  AB  aircraft  body  type 

♦Aircraft  arc  classified  as  WB  (widebody),  NB  (narrow  body), 
♦GM  (ground  manuvcrablc),  or  TT  (tactical). 

/  WB,  NB,  GM,  TT/ 

ENR(AF)  enroute  airfields 
SINCLUDE  cnroulc.txt 

/; 
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♦moving  cargo  requirements  for  units  given  origin  and  destination 
♦origins  and  destinations  are  not  centroids 
PARAMETER  TONSUE(UNIT.I.K) 

SINCLUDE  movingl.txt 

♦moving  passenger  requirements  for  given  origin  and  destination 
♦origins  and  destinations  are  not  centroids 
PARAMETER  TPAX(UNIT,I,K) 

SINCLUDE  movcpal  1  txt 

♦iiiog  capacit)'  for  all  airfield  before  aggregation 
TABLE  MOG(AF,AB) 

SINCLUDE  mogcap.txt 

scalar  totton  total  cargo  after  aggregation 
tottpax  total  passenger  after  aggregation; 
totton=sum((unit.i,k),tonsue(unit,i,k)); 
tottpax=sum((unit,i,k),tpax(unitJ,k)); 

♦aggregation  parameter 

♦aggregate  cargo  and  passenger  due  to  its  origin  and  centroid 
PARAMETER  REPORT2(UNITJ,L); 

PARAMETER  REPORT22(J,L); 

REPORT2(UNIT.J.L)=SUM((LK)$((XORlG.L(I,J)  EQ  1.0)  AND  (XDEST.L(K,L)  EQ  1.0)), 
TONSUE(UNlT,LK)); 

REPORT22(J,L)=SUM(UNIT,REPORT2(UNIT,J.L)); 

PARAMETER  REPORT3(UNIT,J,L); 

REPORT3(UNIT,J,L)=SUM((I,K)$((XORIG.L(l,J)  EQ  1.0)  AND  (XDEST.L(KX)  EQ  1.0)), 
TPAX(UN1T,I,K)); 

PARAMETER  REPORT6(J,AB); 

PARAMETER  REPORT61(L,AB); 

♦aggregation  of  mog  values  into  centroids  first  origin 
REPORT6(J.AB)=SUM(IS(XORIG.L(LJ)  EQ  1.0),MOG(1,AB)); 

♦now  destinations 

REPORT61(L,AB)=SUM(KS(XDEST.L(K,L)  EQ  1.0),MOG(K,AB)); 

PARAMETER  REPORT4(J,L); 

LOOP(UNlT, 

REPORT4(J,L)S((REPORT2(UNIT,J,L)  GT  0.0)  OR  (REPORT3(UNIT.J.L)  GT  0  0))=1 

); 

PARAMETER  REPORT5(*), 

REPORT5(J)S((SUM(L,REPORT4(J,L)))  GT  0.0  )=1.0; 
REPORT5(L)S((SUM(J,REPORT4(J,L)))  GT  0.0)  =1.0; 

FILE  FI  /movement. txt/; 

FILE  F2  /distance. txt/; 

FILE  F3  /pairs  txt/; 

FILE  F4  /mogcap.dat/; 

FILE  F5  /origin  1  l.txt/; 

FILE  F6  /dcstin  1 1  txt/; 
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FILE  F7  /dislancc.dat/; 
FILE  F8  /routclog.txty; 


PUT  FI, 

OPTION  REPORT2;2:():L 
OPTION  REPORT.'?:2 :():!; 

PUT  "MOVEMENT  REQUIREMENTS"; 

LOOP{(UNIT,J.L)$((REPORT2(UNIT,J,L)  GT  0.0)  OR  REPORTS  (UNIT,  J,L)  GTO.O), 
PUT  /; 

PUT  UNIT.TL:4; 

PUT  J.TL:4; 

PUT  L.TL:4; 

PUT " 

PUT  REPORT2(UNIT,J,L); 

PUT " 

PUT  REPORT3(UNIT,J,L); 

); 


PUT  F2; 

PUT  "DISTANCES  BETWEEN  AIRFIELDS"; 

LOOP  ((AF,AFP)$((ENR(AF)  OR  I(AF)  OR  K(AF))  AND 
(ENR(AFP)OR  I(AFP)OR  K(AFP))) , 

PUT  /; 

PUT  AF.TL:4; 

PUT  AFP.TL:4; 

PUT " 

PUT  DIS(AF.AFP); 

): 


PUT  F3; 

PUT  "origin  destination  pairs  for  route  selection  model"; 
LOOP((J,L)$(REPORT4(J,L)  EQ  1.0), 

PUT  /; 

PUT  J.TL:4; 

PUT  L.TL:4; 

); 

PUT  F4; 

PUT"  WB  NB  GM  TT"; 
LOOP(ENR$(REPORT5(ENR)  NE  1.0), 

PUT  /; 

PUT  ENR.TL:4; 

PUT "  "; 

PUT  MOG(ENR,*WB'):2:0; 

PUT"  "; 

PUT  MOG(ENR.'NB’):2:0; 

PUT"  "; 

PUT  MOG(ENR,’GM');2;0; 

PUT"  "; 

PUT  MOG(ENR.'TT'):2:0; 

): 


59 


L00P(J$(REP0RT5(J)  EQ  10), 

PUT/; 

PUT  J.TL:4, 

PUT  ” 

PUT  REPORT6(J,'WB’):2:0; 

PUT" 

PUT  REPORT6(J,'NB'):2:(); 

PUT" 

PUT  REPORT6(J,’GM’):2:0; 

PUT" 

PUT  REPORT6(J,Tr):2:(); 

): 

LOOP(L$(REPORT5(L)  EQ  1.0), 

PUT  /, 

PUTL.TL:4; 

PUT " 

PUT  REPORT61(L,’WB’);2:0: 

PUT" 

PUT  REPORT61(L;NB'):2:0; 

PUT" 

PUT  REPORT61(L,'GM'):2:0; 

PUT  " 

PUT  REPORT61(L,TT'):2:0, 

); 

PUT  F5; 

PUT  "ORIGIN  CENTROID  AIRFIELDS"; 
LOOP(J$(YORIG.L(J)  EQ  1.0), 

PUT  /; 

PUT  J  TL:4; 

): 


PUT  F6; 

PUT  "DESTINATION  AIRFIELDS"; 
LOOP(L$(YDEST.L(L)  EQ  !.()), 

PUT  /; 

PUT  L.TL:4; 

), 

PUT  F7; 

PUT  "*DISTANCES  BETWEEN  AIRFIELDS"; 

PUT  /: 

PUT 

LOOP  ((AF,AFP)S((ENR(AF)  OR  1(AF)  OR  K(AF))  AND 
(ENR(AFP)  OR  I(AFP)  OR  K(AFP)))  , 

PUT  /; 

PUT  AFTL:4; 

PUT"  . 

PUT  AFP.TL;4; 

PUT"  = 
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PUT  DIS(AF,AFP): 

): 

PUT  f  : 

PUT  7:": 

PUT  F8; 

F8.AP-1; 

PUT  “CENTROID  MODEL  RESULTS"; 

PUT  /; 

PUT  “DATE 

PUT  SYSTEM.DATE; 

PUT/; 

PUT  “TIME  : 

PUT  SYSTEM.TIME; 

PUT  /; 

PUT  “ORIGIN  CENTROID  AIRFIELDS”; 

PUT  /; 

LOOP(J$(YORIG.L(J)  EQ  1.0), 

PUT  J.TL;4; 

PUT  /; 

); 

PUT  "DESTINATION  CENTROID  AIRFIELDS"; 

PUT  /; 


LOOP(L$(YDEST.L(L)  EQ  1.0), 

PUT  L.TL:4; 

PUT  /; 

): 

PUT  "MOG  VALUES  FOR  THE  AIRFIELDS"; 
PUT  /; 

PUT"  WB  NB  GM  TT“; 
LOOP(ENR$(REPORT5(ENR)  NE  1.0), 

PUT  /; 

PUT  ENR.TL:4; 

PUT  “  "; 

PUT  MOG(ENR,'WB');2:0; 

PUT"  "; 

PUT  MOG(ENR,'NB’):2:0; 

PUT"  "; 

PUT  MOG(ENR;GM’):2:0; 

PUT"  "; 

PUT  MOG(ENR,TT’):2:0; 

); 

LOOP(J$(REPORT5(J)  EQ  1.0), 

PUT  /; 

PUT  J.TL;4: 
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PUT  ” 

PUT  REPORT6(j;WB'):2:(); 

PUT” 

PUT  REP0RT6(J;NB’):2;(); 

PUT” 

PUT  REP0RT6{J;GM‘):2:(): 

PUT”  ”; 

PUT  REPORT6(j;Tr);2:(): 

): 

L00P(L$(REP0RT5(L)  EQ  l  .O), 

PUT  /; 

PUT  L.TL:4; 

PUT  ” 

PUT  REP0RT61(L;WB’):2:(); 

PUT”  ”; 

PUT  REPORT61(L;NB’):2:0; 

PUT”  ”: 

PUT  REP0RT6UL;GM’):2:(): 

PUT” 

PUT  REPORT61(L,Tr):2:0; 

): 

DISPLAY  0RIG.L,XDEST.L,T0NSUE,TPAX,REP0RT2,REP0RT3,REP0RT4,REP0RT6,REP0RT61; 
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APPENDIX  C.  SENSITIVITY  ANALYSIS  FOR  ROUTES 


STITLE  COLUMN  GENERATION 
SSTITLE  Yasin  TURKER  29  APR  1995 

* - GAMS  AND  DOLLAR  CONTROL  OPTIONS - 

SOFFUPPER  OFFSYMLIST  OFFSYMXREF 

OPTIONS 

LIMCOL  =  0  ,  LIMROW  =  0  ,  SOLPRINT  =  OFF  ,  DECIMALS  =  2 
RESLIM  =  18000,  ITERLIM  =  400000,  OPTCR  =0.1,  SEED  =  3141 


*read  candidate  routes 
SETS 

SINCLUDE  gamscol.dat 

DELR(GEN)  additional  routes  for  column  generation 
DDR(GEN) 

DRR(GEN) 

DRI(GEN,AF) 

DRK(GEN,AF) 

DRE(GEN,AF) 

DADJACENT(GEN,AF,AF) 

DRA(GEN,A); 


DELR(GEN)$(CAN(GEN))=YES$(NOTR(GEN)); 

DDR(GEN)$(DRCAN(GEN))=YES$(NOTDR(GEN)); 

DRR{GEN)$(RRCAN(GEN))=YES$(NOTRR(GEN)); 

DRI(GEN.AF)$(RlCAN(GEN,AF))=YES$(NOTRI(GEN,AF)); 

DRK(GEN,AF)$(RKCAN(GEN,AF))=YES$(NOTRK(GEN,AF)); 

DRE(GEN,AF)S(RECAN(GEN,AF))=YES$(NOTRE(GEN,AF)); 

DADJACENT(GEN,AF,AFP)$(ADJCAN(GEN,AF,AFP))=YES$(NOT  ADJACENT(GEN,AF,AFP)); 
DRA(GEN,A)$(RACAN(GEN,A))=YES$(NOTRA(GEN,A)); 

MAXLEG(GEN)$DELR(GEN)=SMAX((AF,AFP)$DADJACENT(GEN,AF,AFP),DIST(AF,AFP)); 

LOOP((GEN,A)$(DELR(GEN)  *  DRA(GEN,A)), 

CURRENT(AF)=NO; 


*  current  is  set  to  the  first  airfield  on  the  route 

CURRENT(AF)=YES$((DDR(GEN)  *  DRI(GEN,AF))  +  (DRR(GEN)  *  DRK(GEN,AF)) ); 
NEXT(AF)=NO; 

NEXT(AF)=SUM(CURRENT,DADJACENT(GEN,CURRENT,AF)); 

FLTTIME(A,GEN)=SUM((CURRENT,NEXT),DIST(CURRENT,NEXT))/(100*SPD(A)); 

WI(AF)=CURRENT(AF)*I(AF); 

WK(AF)=CURRENT(AF)*K(AF); 

CTIME(A.GEN,NEXT)=GTIME(A,’ONLD')$(CARD(WI)  GT  0)  + 
GTIME(A,*OFFLD’)$(CARD(WK)  GT  0)  + 
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SUM(CURRENT,D1ST(CURRENT,NEXT))/SPD(A); 

KEND=1; 

LOOP(AFP$KEND, 

CURRENT(AF)=NO; 

CURRENT(AF)=YES$NEXT(AF). 

NEXT(AF)=NO; 

NEXT(AF)=SUM(CURRENT,DADJACENT(GEN,CURRENT,AF)); 

FLTT1ME(A,GEN)=FLTTIME(A,GEN)+ 

SUM((CURRENT,NEXT),DIST(CURRENT,NEXT))/(100*SPD(A)); 

CT1ME(A,GEN,NEXT)=SUM(CURRENT,CTIME(A,GEN,CURRENT))+ 

GTIME(A,'ENR’)+ 

SUM(CURRENT,DIST(CURRENT,NEXT))/SPD(A); 

*  quit  inner  loop  if  NEXT  is  empty 
KEND=  CARD(NEXT); 

); 


MAXLOAD(A.GEN)  $  (DELR(GEN)  *  DDR(GEN)  *  DRA(GEN,A))  = 
ACLOAD(A,'r)$(MAXLEG(GEN)  LE  ACRANGE(A,T'))  + 

( ACLOADC  a;  1  ’)+SLOPE(  A,'2')*(  ACRANGE(A,'  1  ')-MAXLEG(GEN))) 
$(MAXLEG(GEN)  GE  ACRANGE(A,T’)  AND 
MAXLEG(GEN)  LE  ACRANGE(A,'2'))  + 
(ACL0AD(A.'2')+SL0PE(A,'3')*(ACRANGE(A,’2')-MAXLEG(GEN))) 
$(MAXLEG(GEN)  GE  ACRANGE(A,'2')  AND 
M AXLEG(  GEN)  LE  ACRANGE( A.G '))  + 

(ACLOAD(A,'3')+SLOPE(A,'4')*(ACRANGE(A,'3')-MAXLEG(GEN))) 
$(MAXLEG(GEN)  GE  ACRANGE(A,'3’)  AND 
MAXLEG(GEN)  LE  ACRANGE(A,'4'))  + 
(ACLOAD(A,'4’)+SLOPE(A,'5’)*(ACRANGE(A.'4')-MAXLEG(GEN))) 
$(MAXLEG(GEN)  GE  ACRANGE(A,'4')  AND 
MAXLEG(GEN)  LE  ACRANGE(A,'5'))  + 
(ACL0AD(A.'5')+SL0PE(A,'6’)*(ACRANGE(A;5')-MAXLEG(GEN))) 
$(MAXLEG(GEN)  GE  ACRANGE(A,'5')  AND 
MAXLEG(GEN)  LE  ACRANGE(A,'6'))  + 
(ACLOAD(A.’6')+SLOPE(A.'7')*(ACRANGE(A,'6')-MAXLEG(GEN))) 
$(MAXLEG(GEN)  GE  ACRANGE(A,'6')  AND 
MAXLEG(GEN)  LE  ACRANGE(A,'7'))  ; 

SONMIXED 


*  first  deliver*'  routes 

KCONSUME(A,GEN.TP,T)$(  (ORD(TP)  LE  ORD(T))  AND  DELR(GEN)  AND 
(DDR(GEN)*DRA(GEN,A))  )=SUM(K$DRK(GEN,K),  CTIME(A,GEN,K)/24); 

KCONSUME(A,GEN,TPT)$(  (ORD(TP)  LE  ORD(T))  AND  DELR(GEN)  AND 
(DDR(GEN)*DRA(GEN,A))  )=MIN(  ORD(T)-ORD(TP)+l  ,  KCONSUME(A,GEN,TP,T)); 

*  now  recovery-  routes 

KCONSUME(A,GEN,TP,T)$(  (ORD(TP)  LE  ORD(T))  AND  DELR(GEN)  AND 
(DRR(GEN)*DRA(GEN,A))  )=SUM(1$DR1(GEN,I),  CTIME(A,GEN,I)/24); 


64 


KCONSUME(A.GEN.TP.T)$(  (ORD(TP)  LE  ORD(T))  AND  DELR(GEN)  AND 
(DRR(GEN)*DRA(GEN.A)) )  =MIN(  ORD(T)-ORD(TP)+ 1  ,  KCONSUME(A,GEN,TP,T)); 

♦+++++++++++++++++++  dynamic  set  DECLARATION  ++++++++++++++++++ 

*  DSETX(U,A,GEN,T)  is  used  to  control  the  allowable  combination 

*  of  unit,  aircraft,  route,  and  start  time  for  the  following 

*  decision  variables:  X(U,A,GEN,T),  TONSUE(U,A,GEN,T)  and 

*  TPAX(U.A.GEN.T).  A  valid  combination  implies  that  the  following 

*  must  be  satisfied: 

*  ( 1 )  it's  a  dcliver>'  route  (2)  aircraft  A  can  fly  route  R 

*  (3)  start  time  >=  unit's  ALD  (4)  delivery  date  <=  RDD  +  MAXLATE 

*  (5)  positive  movement  requirement  (6)  Maximum  aircraft  payload  for 

*  the  route  must  be  greater  than  25%  its  full  load  capability;  this  is 

*  for  cost  efficiency.  (7)  Suppply  of  aircraft  type  >=  0 

DSETX(U,A,GEN,T)$DELR(GEN)  =  YES  $(  1$DDR(GEN)  AND  1$DRA(GEN,A)  AND 
(ORD(T)  GE  ALD(U))  AND 

(ORD(T)  LE  (RDD(U)  -  ROUND(SUM(K$DRK(GEN,K),CTIME(A,GEN,K))/24) 

+  MAXLATE))  AND 

( l$SUM((I.K)$(DRl(GEN,l)*DRK(GEN,K)),MOVEUE(U,I,K))  OR 
1  $SUM((I.K)$(DRI(GEN,I)*DRK(GEN,K)), MOVER AX(U,I,K)) ) 

AND  {l$(MAXLOAD(A.GEN)  GT  (0.25  *  ACLOAD(A.'T")) )) 

AND  SUM  (TP$(ORD(TP)  LE  ORD(T)),  SUPPLY(A,TP)) ) ; 


*  DSETDY(A,GEN,T)  is  used  to  control  the  allowable  combination 

*  of  aircraft,  route,  and  start  time  for  the  decision  variable 

*  Y(A,GEN,T).  An  allowable  combination  is  one  in  which  the 

*  aircraft-route  is  compatible  and  that  the  supply  of  aircraft 

*  type  to  date  is  greater  than  0, 

DSETY(A,GENT)$DELR(GEN)  =  YES  $(  1$DRR(GEN)  AND  1$DRA(GEN,A)  AND 
SUM  (TP$(ORD(TP)  LE  ORD(T)),  SUPPLY(A,TP)) ); 

*definc  priceout  parameters 


PARAMETER  TOTX(U,A,GEN,T)  price  out  value  for  x; 
TOTX(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=0.0; 

"■computations  for  TOTX 

TOTX(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=TOTX(U,A,GEN,T)-h 
SUM(1$(DSETX(U,A,GEN,T)"‘DRI(GEN,I)),ACBALI.M(A,I,T)); 

LOOP(TP, 

TOTX(U.A.GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=TOTX(U,A,GEN,T)- 
SUM(K$(1$DSETX(U,A,GEN,T)  AND  1$DRK(GEN,K)  AND 
((ORD(T)  +  ROUND(CTIME(A.GEN,K)/24))  EQ  ORD(TP))), 

ACBALK.M(A,K,TP)) ); 

TOTX(U.A.GEN.T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=T0TX(UAGEN,T)- 
(FLTT1ME(A,GEN)*ACURATE.M(A)); 

LOOPfTP. 

TOTX(U.A.GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN)  AND  (ORD(T)  LE  ORD(TP))) 
TOTX(U,A,GEN,T)+(KCONSUME(A,GEN,T,TP)*ACCONSUME.M(A,TP)) ); 
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LOOP(TP. 

TOTX(U.A.GEN.T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=TOTX(U,A,GEN,T)+ 
SUM(AF$(DRJ(GEN,AF)»DDR(GEN)), 

MOGREQ(AF,A)*(GTIME(A,'ONLD')/24)*MOGUTILITY.M(AF,T))+ 
SUM(AF$(1$DRK(GEN,AF)  AND  1$DDR(GEN)  AND 
(ORD(T)  +  ROUND(CTIME(A,GEN,AF)/24)  EQ  ORD(TP))), 
MOGREQ(AF,A)*(GTIME(A,'OFFLD')/24)*MOGUTILITY.M(AF,TP))+ 
SUM(AF$(1$DRE(GEN,AF)  AND  1$DDR(GEN)  AND 
(ORD(T)  +  ROUND(CTIME(A.GEN.AF)/24)  EQ  ORD(TP))), 
MOGREQ(AF,A)*(GTIME(A,'ENR’)/24)*MOGUTILITY.M(AF,TP)) ); 
*computation  for  TONSUE 


PARAMETER  TOTTON(U,A,GEN,T); 

TOTTON(U,A,GEN,T)$(DSETX(U,A,GENT)  AND  DELR(GEN))=0.0; 
TOTTON(U,A,GEN,T)$(DSETX(U,A,GEN.T)  AND  DELR(GEN))=TOTTON(U,A,GENJ)+ 
SUM((I.K)$(ACSQFT(A)  AND  1$DSETX(U,A,GEN,T)  AND  1$DRI(GEN,1)  AND 
1$DRK(GEN,K)  AND  l$MOVEUE(U,I,K)),REQMTUE.M(U,I,K)); 

TOTTON(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=TOTTON(U,A,GEN,T)+ 
SUM((I,K)$(l$MOVEUE(U,I,K)  AND  ACCARGO(A.'OUT')  AND 
1$DSETX(U,A,GEN,T)  AND  1$DRI(GEN,I)  AND  1$DRK(GEN,K)), 
OUTREQMT.M(U.I,K)); 

TOTTON(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=TOTTON(U,A,GEN,T)+ 
SUM((l,K)$(l$MOVEUE(U,I,K)  AND  ACCARGO(A,'OVER')  AND 
1$DSETX(U,A,GEN.T)  AND  1$DRI(GEN,I)  AND  1  $DRK(GEN,K)), 
OVERREQMT.M(U,l,K)), 


^includes  cost  coefficient 
LOOP(K, 

TOTTON(U.A.GEN,T)$(ACSQFT(A)  AND  1$DSETX(U.A.GEN,T)  AND  1$DELR(GEN)  AND 
l$DRK(GEN,K))=TOTTON(U,A,GEN,T)- 

LATEPEN(U,'UE’)*MAX((),(ORD(T)+ROUND(CTIME(A,GEN,K)/24)-RDD(U))) ); 


*computations  for  variable 
PARAMETER  TOTY(A,GEN,T); 

TOTY(A,GEN,T)$(DSETY(A,GEN,T)  AND  DELR(GEN))=().0; 

LOOP(TP, 

TOTY(A,GEN,T)$(DSETY(A,GEN,T)  AND  DELR(GEN))=TOTY(A,GEN,T)- 
SUM(I$(1$DR1(GEN,1)  AND 

((ORD(T)+ROUND(CTIME(A,GEN,l)/24))  EQ  ORD(TP))),ACBALI.M(A,I,TP)) ); 

TOTY(A,GEN,T)$(DSETY(A,GEN,T)  AND  DELR(GEN))=TOTY(A,GEN,T)+ 
SUM(K$(DSETY(A,GEN,T)*DRK(GEN,K)),ACBALK.M(A,K,T)); 

TOTY(A,GEN,T)$(DSETY(A,GEN,T)  AND  DELR(GEN))=TOTY(A,GEN,T)+ 
(FLTTIME(A,GEN)»ACURATE.M(A)); 

LOOP(TP, 

TOTY(A,GEN,T)$(DSETY(A,GEN,T)  AND  DELR(GEN)  AND  (ORD(T)  LE 
0RD(TP)))=T0TY(A,GEN,T)+ 

(KC0NSUME(A,GEN,TTP)*ACC0NSUME.M(A,TP))  ); 
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LOOP(TP. 

TOTY(A.GEN,T)$(DSETY(A,GEN,T)  AND  DELR(GEN))=TOTY(A,GEN,T)+ 

SUM(AF$(1$DSETY(A.GEN,T)  AND  1$DRE(GEN,AF)  AND  1$DRR(GEN)  AND 
l$(ORD(T)  +  ROUND(CTIME(A,GEN,AF)/24)  EQ  ORD(TP))), 
M0GREQ(AF,A)*(GT1ME(A,'ENR')/24)*M0GUT1LITY.M(AF,TP))); 


♦computations  for  variable  TPAX 
PARAMETER  TOTTPAX(U,A,GEN,T); 

TOTTPAX(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=0.0; 
TOTTPAX(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=TOTTPAX(U,A,GEN,T)+ 
SUM((I,K)$(1$DSETX(U,A,GEN,T)  and  1$DRI(GEN,I)  and  1$DRK(GEN,K)AND 
l$MOVEPAX(UJ,K)),REQMTPAX.M(U,I,K)); 

♦include  cost 
LOOP(K, 

TOTTPAX(U,A,GEN,T)$(DSETX(U,A,GEN,T)*DELR(GEN)  *  DRK(GEN,K))= 
TOTTPAX(U,A,GEN,T)- 

(LATEPEN(U.'PAX’)*MAX(0,ORD(T)+ROUND(CTIME(A,GEN,K)/24)-RDD(U))) ); 


♦warmup  starting 

X. M(U,A,GEN.T)$(DSETX(U,A,GEN.T)  AND  DELR(GEN))=  1 .0; 

Y. M(A.GEN,T)$(DSETY(A,GEN.T)  AND  DELR(GEN))=1.0; 
TONSUE.M(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=1.0; 
TPAX.M(U,A,GEN.T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=1.0; 
ACWE1GHT.M(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=1.0; 
ACSPACE.M(U.A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=1.0; 
ACPAX.M(U,A,GEN,T)$(DSETX(U,A,GEN,T)  AND  DELR(GEN))=1.0; 

♦define  new  route  set  base  to  price  out  values 
LOOP((U.A,T), 

R(GEN)$(DSETX(U,A,GEN.T)  *  DELR(GEN))=YES$((TOTX(U.A,GEN,T)  LT  0.0)  OR 
(TOTTPAX(U.A.GEN,T)  LT  0.0)  OR  (TOTTON(U,A,GEN,T)  LT  0.0)  OR  (TOTY(A,GEN,T)  LT  0.0)) ); 
DR(GEN)$DDR(GEN)=YES$R(GEN); 

RR(GEN)$DRR(GEN)=YES$R(GEN); 

RI(GEN,AF)$DRI(GEN,AF)=YES$R(GEN); 

RK(GEN,AF)$DRK(GEN,AF)=YES$R(GEN); 

RE(GEN,AF)$DRE(GEN,AF)=YES$R(GEN); 

ADJACENT(GEN,AF,AFP)$DADJACENT(GEN,AF,AFP)=YES$R(GEN); 

RA(GEN,A)$DRA(GEN,A)=YES$R(GEN); 

RESULT('NEWR’)=CARD(R); 

DISPLAY  RESULT; 


SOLVE  AIRLIFT  USING  RMIP  MINIMIZING  Z; 
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APPENDIX  D.  FLOW  DIAGRAM  OF  ROUTE  SELECTION  MODEL 


START 


TPI  ni).  TXT  (1) 
()RIC;iN.TXT  (2) 
DICSTINA.TXT  (2) 


AGGREG.  PAS  (3) 

TPFDD.  TXT 
ORIGIN.TXT 
DESTINA.TXT 


#c)K(:i:ntroids 

(Kl) 


DESTSET  GMS 

DES'iiNAi  rx^r 

n  c>i  c:i:n  i  Roins  (K) 


!vio<;c:ap.  txt 

KNROin  i:.  I X I 


UNITJ.  TXT 
MOVINGl.  TXT 
MOVEPAll  TXT 
ORIGIN]  TXT 
DESTINAl  TXT 
UESQFTl  TXT 
C ARGOT  I  TXl^ 
ALDI  TXT 
RDDI  TXT 
ROUTELOG  TXT  (9) 

AGGREG  GMS  (4) 

UNITL  DAT 
MOVINGl  DAT 
MOVEPAll  DAT 
MOGCAP.  TXT 
ENROIJTE.  TXT 

MOVEMENT  TXT 
PAIRS  TXT 
DISTANCE.  TXT 
MOGCAP  DAT 
ORIGIN  TXT 
DESTINA.  TXT 
DISTANCE.DAT 
ROUTELOG.  TXT 


AIRFIELD.  SET 
DISTl.  TXT 
#OFCENTROIDS  (K) 


ORIGSET.  GMS 

AIRFIELD.  SET 
DISTl.  TXT 
#OFCENTROIDS  (K2) 
ORIGIN!  TXI' 


SECAGGRE.  PAS  (5) 

ALDI  TXT' 

RDDI  TXT 
MOVEMENT  TXT 
CARGOPl  TXI' 
UESQFTl.  TXT 
ORIGIN  TXT 
DESIINA  TXT 

CONT, 
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ENROUTE.  TXT 
CREW.  TXT 
AIRFIELD.  TXT 
AIRCRAFT.  TXT 
LAMBDA.  TXT 


SECAGGRE  PAS 

IJNITNAME.DAT 
MOVEUE.  DAT 
MO  YEP  AX  DAT 
IJESQPT.  DAT 
CARGO P.  DAT 
ALDT.  DAT 
RDDT.  DAT 
ROIJTELOG.  TXT 


ROUTESEL  PAS  (6) 

PAIRS.  TXT 

ENROUTE.  TXT 
I  CREW.TXT 
AIRFIELD.  TXT 
/  AIRCRAFT.  TXT 
DISTANCE.TXT 
LAMBDA.  TXT 

GAMSREG.  DAT 
ROWE.  TXl' 

DAVE.  DAT 
USER.  DAT 
ROIJTELOG.  TXT 


THROUHPUT  GMS  (7) 

GAMSREG.  DAT 
MOVEUE.  DAT 
MOVE  PAX.  DAT 
DISTANCE.DAT 
MOGCAP.  DAT 
UNITNAME.  DAT 
IJESQrr.  DAT 
CARGOP  DAT 
ALD.  DAT 
RDD.  DAT 


COLUMN.  GMS  (8) 
GAMSCOL  DAT 


STOP 
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1  Text  file  separated  with  tab  and  generated  in  EXCEL. 

2  Preselected  origin-destination  airfields  which  will  be  considered  during  the 
aggregation  If  the  user  wants  to  include  all  the  airfields  that  mentioned  in  TPFDD, 
then  enter  “No”  for  the  including  file  question. 

3  Unit  aggregation  based  on  preselected  origin  destination  airfields  ,  minimum 
movement  requirements  and  the  difference  between  the  RDD  values. 

4  Centroid  model  for  origin  and  destination  airfields  written  in  GAMS. 

5.  Unit  aggregation  based  on  centroid  airfields. 

6.  Route  selection  model. 

7  Optimization  model. 

8  Sensitivity  Analysis. 

9  Text  file  contains  all  the  input  and  aggregation  information  based  on  this  run. 
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